The New Heavy Mesons: A Status Report 



Eric S. Swanson* 
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A survey of the experimental, phenomenological, and theoretical status of the new heavy mesons 
is presented. States discussed are the B c , h c , rj c , D s (2317), L> 3 (2460), X(3872), X(3940), Y(3940), 
Z(3930), and y(4260). Quark models for spectra, strong decays, and hadronic interactions are 
reviewed and used to interpret the new states. New results for strong decay models, bound state 
decays, mesonic molecules, properties of the X(3872), and the chiral doublet model are also pre- 
sented. 
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I. INTRODUCTION 



More than 30 years after the November Revolution, charmonium spectroscopy continues to surprise and challenge. 
A new era began in April of 2003 when BaBar announced the discovery of the enigmatic Z) s (2317). Since then CLEO, 
Belle, Fermilab, and BES have joined the scrum, and the number of new states has risen to the double digits. 

The formerly comfortable world of heavy quark spectroscopy is being severely tested by new experiment [l]. This 
report summarises the experimental status of eleven new open and hidden charm mesons 1 , the challenges these make 
to our understanding of strong QCD, and the theoretical attempts to meet these challenges. 

This report begins with a summary of the new eleven and a review of what is expected in charmonium spectroscopy. 
For the latter, several models (new and old) are compared to the data, and three models of strong decays are examined. 
This serves to illustrate what is known (and not known) before examining the new charmonia states. A similar 
introduction to heavy-light systems is made in Sections IV. A and IV. B. 



A. Summary 



Table [I] summarises the contents of this report for the impatient. A selection of interpretations and theoretical ideas 
are listed in the 'comments' column. Theory and more experimental detail can be found in the main text. 



B. Charmonia 



This section describes the current status of charmonium and strong decay (mostly constituent quark and lattice 
gauge theory). This is in an effort to set the stage for the new hadrons to be described in subsequent sections and to 
establish the degree of reliability which we may attach to those expectations. 

The utility of the nonrelativistic limit in describing charmonia has been recognised since the days of the November 
Revolution^, 0- The subsequent advent of the operator productQ, nonrelativistic and potential nonrelativistic 
QCD 6] formalisms has considerably strengthened the connection of heavy quark phenomenology with QCD 2 . For 
example, the hierarchy of scales m c >> m c a s (m c ) >> m c a 2 s {m c ) permits one to choose a cutoff which removes 
transverse gluons from the effective description of a meson. Thus an efficient potential description and Fock space 
expansion of charmonium is feasible |7j. 



1 States not discussed here include the X(1835), putative glueballs, hybrids 7ri(1400), 7ri(1600), and 7r(1800), and baryons. 

2 It also reveals that the charm quark is uncomfortably light. 
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TABLE I: New Heavy Mesons 



state 


mass (MeV) 


width (MeV) 


production/decay mode 


comments 


ref 


h c 


3524.4 ±0.6 ±0.4 


- 


V(25)^^ u / tc -,(77)(7»7c) 


~ CQM / tests spin dependence 


CLEO[87] 


rfc 


3654 ± 6 ± 8 


< 55 


B — > Kr\ c — ► KK s Kn 


~ CQM / tests hyperfine splitting 


Belle [90] 




3642.9 ± 3.1 ± 1.5 


6.3 ± 12.4 ±4.0 


e + e~ —* ri' c J/ip 




CLEO[93] 




3630.8 ± 3.4 ± 1.0 


17.0 ±8.0 ±2.5 


77 — > — ► K$K~k 




BaBarf94] 


X(3872) 


3872.0 ±0.6 ±0.5 


< 2.3 95% C.L. 


B -> O -f K-k-kJ/iP 


molecule, cusp, tetraquark 


Belle [32J 




3873.4 ± 1.4 


- 


B ^ KX ^ K-k-kJ/iP 




BaBarl35] 




- 


- 


B — > X — > nn-Kj/tp 




Belle [43] 




- 


- 


B -» X -> 7J/^ 




Belle [43J 




3871.3 ±0.7 ±0.4 


— 


pp — > X — > TTnJ/lp 




CDF [331 




3871.8 ± 3.1 ± 3.0 


— 


pp — > X — > TrnJ/ip 




D0[34l 




avg = 3871.9 ±0.5 










X(3940) 


3943 ± 6 ± 6 


< 52 


e + e~ — ► J/ipX — ► J/ipDD* 


Xcii ??" / needs confirmation 


Belle [96] 


F(3940) 


3943 ± 11 ± 13 


87 ± 22 ± 26 


B — > — > KirnirJ/ij) 


needs confirmation 


Bellcll05] 


Z(3930) 


3931 ± 4 ± 2 


20 ± 8 ± 3 


77 — > Z — > D£> 


Xc2 / » CQM 


Bellell07] 


Y (4260) 


4259 ± 8 ± 4 


O O 1 C\ O 1 f 

88 ± 23 ± 5 


e T e — > 7/sfl^ — * -yisnJ/W^ 


hybrid:/ needs conhrmation 


BaBarll09J 


DJ2317) 


2317.3 ± 0.4 ± 0.8 


< 10 


e + e~ — > D s (2317) — » D s 7r u 


molecule, tetraquark, shifted cs 


BaBarll21] 




2319.8 ±2.1 ±2.0 


w 


B -> L>D S (2317) -> _D£> s7 r° 




Bellell39] 




2318.5 ± 1.2 ± 1.1 




L> s (2317) -> L> s tt 




CLE011321 


D s (2460) 


2463.6 ± 1.7 ± 1.2 


< 7 90% C.L. 


D s (2460) -> D>° 


molecule, tetraquark, shifted cs 


CLEQ11321 




2458.0 ± 1.0 ± 1.0 


resolution 


A, (2460) -> D s 7r°7 




BaBarll34l 




2459.2 ± 1.6 ±2.0 


w 


S -> DA,(2460) -> DD* a Tv°, DD S ^ 




Bellell39l 


L> s (2630) 


2632.6 ± 1.6 


< 17 90% C.L. 


D s -» and D s?? 


artefact 


SELEX1168] 


B c 


6285.7 ±5.3 ±1.2 


0.474 ±0.07 ±0.33 ps 


pp^ B c —> J/f/'7r ± 


« CQM 


CDF1180] 



1. Spin-dependent Interactions 



Typical nonrelativistic quark models contain the following ingredients: nonrelativistic quark kinetic energy, a central 
confining potential, and a variety of spin-dependent (v/c) 2 corrections. Thus one has 



H C QM = K(q) ± K(q) + V conf + V S D (1) 

where if is a quark kinetic energy, V con f is the central confining potential, and Vsd is the spin-dependent potential. 

The form of the spin-dependent interaction has been computed by Eichten and FeinbergQ at the tree level using 
the Wilson loop methodology. This result was extended by Pantaleone, Tye, and NgQ in a one-loop computation. 
The net result is: 



V SD {r) 



( a i 

\4m 2 



4m? 



1 de 
r dr 



± 



Ylm q mq 



( 3cr 9 



r dr 

(T q ■ (Tg 
(T q — (Tq 

m„mn 



V 3 



2m q mq 



Ylm q rriq 



ldV2 

r dr 



CTqVi 



L 



(2) 



Here e = e(r) is the static potential, L = L 9 = — L q , r = \r\ = |r g — r<j | is the QQ separation and the Vi = Vi(m 
are determined by electric and magnetic field insertions on quark lines in the Wilson loop matrix element. V5 vanishes 
in the equal quark mass case. As shown by Gromes [lfj covariance under Lorentz transformation leads to a constraint 
between the SD potentials, 



e'{r) = Vi{r)-V[{r), (3) 

In a model approach, QQ interactions are typically derived from a nonrelativistic reduction of a relativistic current- 
current interaction 
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H int = \J d 3 xd 3 yJ r (x)K{x - y)J r {y) (4) 

where Jr = ipTip specifies a quark current. 

As far as long range potentials are concerned only time-like vector or scalar currents are relevant since other currents 
do not yield static central potentials|l2j. Vector confinement interactions are motivated as natural extensions of one- 
gluon-exchange. More fundamentally, the instantaneous nonperturbative interaction of QCD in Coulomb gauge has 
the Lorentz structure V° <S>V°. Scalar interactions are also often considered when model building for reasons to be 
discussed shortly. 

Nonrelativistic reductions of vector or scalar interactions yield: 



Vi = 0, V 2 = e, V 3 = e'/r - e" , V 4 = 2V 2 e vector 

Vi = -e, V 2 = V 3 = Vi=0 scalar. (5) 

Determining the Dirac structure of the long range potential yields important information on the nonperturbative 
gluodynamics which gives rise to confinement. This can be accomplished by examining spin-dependent mass splittings 
in quarkonia, in particular the ratio of P-wave splittings defined by 



1++ - 0++ (6) 



is a useful diagnostic |Tl). 

Perturbative P-wave quarkonium splittings are given by 



4 

where 



2++ = E Q -^T + L + - A S (7) 
5 4 

1++ = E Q + 2T - L+ -S (8) 
4 

0++ = E Q -4T-2L+^S (9) 

1 + - = E -^-S (10) 



12m 



' <^3) (11) 



a 2 
1 



L = —{ T + ^ + -^) (12) 
Tfi Zr r r 



3m 2 

Using the results given in Eq. [S] gives vector and scalar spin splittings ratios of 



(V*). (13) 



and 



'16 o, , 14 b\ 



p(V) = ^ (14) 



/16 ct a _ b\ 

p(s) = )": 3 _ (15) 

The measured values for this parameter are p(xc) — 0.48(1) and p(xb) = 0.66(2). Notice that the results given here 
imply that p(QQ) — > 4/5 in the heavy quark limit. Thus, although the motion is in the right direction, even the \b 
system is not yet heavy by this metric. 



5 



Fig. n displays the splitting ratio versus the RMS meson radius along with the experimentally measured ratios for 
the ip and T systems. It is clear that the vector confinement hypothesis cannot explain the data. Alternatively, the 
scalar hypothesis decreases with increasing RMS radius as required. Arrows in the figure show RMS radii for the Xc 
and Xb systems computed in a simple nonrelativistic model, indicating that scalar confinement is also quantitatively 
reliable. 

1.2 
1 

0.8 



a. 0.6 



0.4 
0.2 


1 2 3 4 5 

<r> (GeV" 1 ) 

FIG. 1: Mass Splitting Ratio vs. rjjMS- Experimental values are given as horizontal lines. 

The hyperfine interaction V4 may be probed by comparing the mass splitting between the P-wave system and the 
axial vector quarkonium state. Specifically, perturbative tensor and spin-orbit splittings are not present in the 'centre 
of gravity' of the P-wave system: 

M cog = ~ (M(0++) + 3M(1++) + 5M(2++)) = E Q + ]s. (16) 
y 4 

Thus 

M cog - M (1+-) = S. (17) 

For one gluon exchange S oc (S(r)) and since this expectation value is zero in P-waves, very little axial-cog mass 
splitting is expected (little because the delta function is smeared by zitterbewegung effects). Alternatively, for a 
vector confining potential S = ancl a nonzero splitting can be expected. The recent discovery of the h c 

permits this splitting to be measured. The results are discussed in Section IlII Al 

Simple spectroscopy teaches us about the v 2 /c 2 structure of the confinement potential, and therefore something 
about the nonperturbative gluodynamics associated with flux tubes. This is nontrivial and worthy of pursuit! 

2. Constituent Quark Model and Lattice Spectra 

The predictions of four representative constituent quark models are given in Tabic ITU These are a simple nonrel- 
ativistic model an updated version of the 'relativised' Godfrey-Isgur model |T3|, a 'relativistic' quark model [Tij. 
and the original Cornell modelQ- 

The simple constituent quark model (labelled BGS in Table [HJl employs a nonrelativistic quark kinetic energy, 
K(q) = m q — V 2 /2m 9 , a 'Coulomb+lincar' central potential, V con f — — 4a s / (3r) + br and a spin-dependent interaction 
that corresponds to the vector one-gluon-exchange + scalar confinement Ansatz in the equal quark mass case: 

V 1 = -br, V 2 = -^, V 3 =4^, V 4 =^5», V 5 = (18) 
3 r r 5 3 
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where S a (r) = (a/y/ir) 3 e a r represents a 'smeared' delta function. Fitting the charmonium spectrum yields: 
a s = 0.5462, b = 0.1425 GeV 2 , m c = 1.4794 GeV, and a = 1.0946 GeV. 

The Godfrey-Isgur model is a "relativised" extension of typical nonrelativistic quark models. In particular, a 

relativistic dispersion relation for the quark kinetic energy K(q) = + V 2 , a QCD-motivatcd running coupling 

a s {r\ & flavor-dependent potential smearing parameter a — > a(m q ,m q ), and additional factors of m q /K(q) are all 
used jl^l Hsj . Just as in the nonrelativistic model, the quark-antiquark potential V qq (p,r) incorporates the Lorentz 
vector one gluon exchange interaction at short distances and a Lorentz scalar linear confining interaction. Note that 
the string tension and quark mass (b = 0.18 GeV 2 and m c = 1.628 GeV) are significantly larger than the values used 
in the nonrelativistic model. Predictions of this model are in the column labelled GI in Table |nj 

Ebert, Faustov, and Galkin[l4| attempt a more formal ('quasipotential') approach to relativising the quark model. 
A nonrelativistic dispersion relation is used but with reduced and total masses determined relativistically (for example, 
the reduced mass is given by [i = E q E q /(E q + E q ) with E q — (M 2 — m? + m 2 )/(2M) and M is the eigenenergy) . 
The interaction employs the complete Coulomb gauge gluon propagator (so that the temporal and spatial vector 
interactions are included), a scalar confinement potential, and a vector confinement potential with an anomalous 
magnetic moment k for the quark. In the nonrelativistic limit, the confinement potential is V con f(r) — br + B where 
the linear potential receives contributions from the vector and scalar confinement terms of (1 — S)br + B and 5br 
respectively where 5 is a mixing coefficient. The coupling is fixed for each quark flavour and is a s = 0.314 for 
charmonia. Additional parameters are determined to be m c — 1.55 GeV, b = 0.18 GeV 2 , B = —0.16 GeV, k = — 1. 
Finally, at 0(v 2 /c 2 ) the spin-dependent interaction reduces to precisely that used in the BGS model discussed above 
(when k = — 1). For example, the long range contribution to V4 is 



V 4 (r) = 4(1 + ac) 2 (1 - 5)- (19) 
r 

which is disfavoured by the data so that that this choice of the anomalous magnetic moment eliminates unwanted 
long range chromomagnetic interactions 0|- The predictions of this model are labelled EFG in Table HTl 

The final model is that of the Cornell group Q, whose original predictions are reproduced in the sixth column of 
Table HT1 The starting point is an assumed nonrelativistic kinetic energy and a quark-quark interaction mediated by 



Hint = \ I ,f\v,f'uo"U)KU -</)/."(//) (20! 

with the kernel given by K(r) — a s /r — |6r and the current given by the colour charge density, p a — ip^T a ip. Of 
course this is a vector confinement model so that the nonrelativistic reduction of Eq. |201 corresponds to 



4 4 a b b 

e(r)=--K(r), V 1= 0, V 2 = --K{r), V 3 = 4^ + -, V A = -8na s S{r) + 4- . (21) 
d 6 r a r r 

Fitted parameters are m c = 1.84 GeV, b = 0.183 GeV 2 , and a s = 0.39. Spin-dependent effects are ignored in the 
original Cornell papers. However, the Cornell group does attempt to include the effects of open charm virtual states 
on the spectrum (this will be discussed in greater detail in Sections ll B 31 and II B 4|l . The resulting masses are shown in 
square brackets in Table [H] We note that the coupled channel approach shifts most of the predicted masses a modest 
amount, except the ifj(3S) and ip{4S) states, for which the naive predictions have been considerably worsened. 

The advent of QCD motivated quark models coincided with the invention of lattice gauge theory and the lattice 
has provided an important counterpoint to models ever since. Recently, lattice computations have become sufficiently 
robust that attempts have been made to compute extensive spectra. Two such computations are presented in the 
last columns of Table ITU These computations have been made in the 'quenched' approximation where the effects of 
virtual quarks are neglected on anisotropic lattices with a lattice volume fixed at (1.6 fm) 3 in Ref. 15] or over lattice 
volumes ranging from (1.5 fm) 3 to (3.3 fm) 3 in Ref. [T^- Continuum extrapolated values are quoted in the table. 

Overall the agreement with the data is quite impressive. Average errors for the models are listed in Table II I II 
Surprisingly, the simple nonrelativistic BGS model is most accurate in describing the spectrum. It is somewhat 
disappointing that lattice has not caught up to models yet, but this may change soon. 

Spectra are not particularly precise determinants of model efficacy. Other observables such as strong and leptonic 
widths and radiative transition rates provide sensitive probes to dynamics and it is important to consider the full set 
of model predictions when making comparisons. 

An example of these considerations concerns the strong and leptonic widths of the ^(3770). The ■0(3770) is generally 
assumed to be the 3 Di cc state; however, mixing with the 2 3 Si Fock component is possible and would modify the 
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TABLE II: Experimental and theoretical spectrum of cc states. States assigned in this review are denoted a . 



State 




BGS[13] GF13] 


1 / i v 1 r-i a l 

hjt (j| t4| 


Cornell lo| 


t~*T> da pori 
G-r--rALjOllo| 


Onenl tb| 


JMi 3 Si) 


3096.87 ± 0.04 


3090 


3098 


3096 


3095 [3095] 


3085 ± 1 


3084 ± 4 




2979.2 ± 1.3 


2982 


2975 


2979 


3095 


3013 ± 1 


3014 ±4 




3685.96 ± 0.09 


3672 


3676 


3686 


3684 [3684] 


3777 ± 40 


3780 ± 43 


^(2 1 S ) 


3637.7 ± 4.4 


3630 


3623 


3588 


3684 


3739 ± 46 


3707 ± 20 


^(3 3 Si) 


4040 ±10 


4072 


4100 


4088 


4110 [4225] 


- 


- 






4043 


4064 


3991 


4110 


- 


- 


V>(4 d S!) 


4415 ± 6 


4406 


4450 


- 


4460 [4625] 


- 


- 






4384 


4425 


- 


4460 


- 


- 


X2(l d P 2 ) 


3556.18 ±0.13 


3556 


3550 


3556 


3522 [3523] 


3503 ± 24 


3488 ±11 


Xi(l 3 Pi) 


3510.51 ±0.12 


3505 


3510 


3510 


3522 [3517] 


3472 ± 9 


3462 ± 15 


Xo(l 3 Po) 

Mi'Pi) 


3415.3 ±0.4 


3424 


3445 


3424 


3522 


3408 


3413 ± 10 


3524 ± l a 


3516 


3517 


3526 


3522 [3519] 


3474 ± 40 


3474 ± 20 


X2(2 J P 2 ) 


3931 ± 5 a 


3972 


3979 


3972 


- 


4030 ± 180 


- 


Xi(2 3 Pi) 




3925 


3953 


3929 


- 


4067 ± 105 


4010 ± 70 


Xo(2 3 Po) 




3852 


3916 


3854 


- 


4008 ± 122 


4080 ± 75 


/i c (2 1 Pi) 




3934 


3956 


3945 


- 


4053 ± 95 


3886 ± 92 


X2(3 d P 2 ) 




4317 


4337 


- 


- 


- 




Xi(3 3 Pi) 




4271 


4317 


- 


- 


- 




Xo(3 3 Po) 
fec^Pi) 




4202 


4292 


- 


- 


- 






4279 


4318 


- 


- 


- 




i/>3(l 3 D 3 ) 




3806 


3849 


3815 


3810 


- 


3822 ± 25 


<M1 3 D 2 ) 




3800 


3838 


3811 


3810 


- 


3704 ± 33 


^(1 3 D!) 


3769.9 ± 2.5 


3785 


3819 


3798 


3810 [3755] 


- 


- 


^(l'Da) 




3799 


3837 


3811 


3810 


- 


3763 ± 22 


^ 3 (2 J D 3 ) 




4167 


4217 


- 


4190 


- 


- 


V 2 (2 3 D 2 ) 




4158 


4208 


- 


4190 


- 


- 


V>(2 3 Di) 


4159 ± 20 


4142 


4194 


- 


4190 [4230] 


- 


- 


?7c 2 (2 1 D 2 ) 




4158 


4208 


- 


4190 


- 


- 


X4(l d F 4 ) 




4021 


4095 


- 


- 


- 


- 


X3(1 3 F 3 ) 




4029 


4097 


- 


- 


- 


4222 ± 140 


X2(1 3 F 2 ) 




4029 


4092 


- 


- 


- 




« C 3(1 Fa) 




4026 


4094 








4224 ± 74 


X4(2 J F 4 ) 




4348 


4425 










X3(2 3 F 3 ) 




4352 


4426 










X 2 (2 3 F 2 ) 




4351 


4422 










fteS^Fa) 




4350 


4424 










^5(1 3 G 5 ) 




4214 


4312 










iMl 3 G 4 ) 




4228 


4320 










V'3(1 3 G 3 ) 




4237 


4323 










r? c4 (l 1 G 4 ) 




4225 


4317 











TABLE III: Average Model Errors. 
BGS GI EFG Cornell CP-PACS Chen" 
0.3% 0.6% 0.4% 1.4% (19% L3%T 



predicted strong and leptonic widths. Allowing this mixing brings a strong decay calculation into agreement with 
the measured width of the ?/;(3770) for a mixing angle of approximately —17°. Assuming that leptonic widths are 
dominated by S-wave wavefunction components then yields the prediction fl3| 



W3770)) 



r e+e -(^(3686)) 

which is in good agreement with the experimental value of 0.12 ± 0.02. 



= 0.10 ± 0.03, (22) 

thy. 
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3. Strong Decay Models and Predictions 



Strong decays of mesons are driven by nonperturbative gluodynamics and are therefore (i) difficult to compute 
(ii) good probes of str ong QCD. Unfortunately, lattice computations of strong decays are very difficult (preliminary 
attempts are in Ref. |2flj) and one is forced to rely on models. The most popular are the l3 Po' modelj2f| with a 
two-body transition operator given by 



and the Cornell model which takes Eq. [2U] seriously. The Cornell model has the advantage of unifying the description 
of the spectrum and decays and completely specifies the strength of the decay. 

Although, extensive comparisons of models have not been made, Ref. [22] finds that 3 Pq amplitudes tend to be 
much larger than those due to perturbative gluon pair production, while Ref. |2^| finds that the 3 P model gives 
better predictions for amplitude ratios than does a model which creates quark pairs with 3 Si quantum numbers. 
Table llVl is an attempt to partially rectify this situation by presenting the predictions of the 3 Po model, the Cornell 
model, and computations of strong widths made in a coupled channel approach to the Cornell model0,0|. The first 
two computations are made in first order perturbation theory, while the latter is nonperturbative. 1 shall refer to 
perturbative Cornell model as the pKp model to highlight this difference. 

The 3 Pq and pKp calculations assume SHO wavefunctions with a universal scale parameter set to (3 = 0.5 GeV; 
meson masses are taken as shown. The 3 Pq model results are those of Ref. 01 with the exception that the decay 
strength has been taken to be 7 = 0.35, giving somewhat better agreement with PDG widths. The pKp computation 
sets K(r) = —3/4br and uses a string tension of b = 0.18 GeV 2 . The last column displays an attempt to extract 
widths from the coupled channel Cornell model (state masses and wavefunctions only approximate those used in the 
other models). Table HVI also lists PDG measured widths and the results of a recent reanalysis of R by SethpEj. 

The pKp and Cornell computations should agree in the event that similar wavefunctions and meson masses were 
employed and the resulting width is not too large. The few cases where it is possible to make a comparison indicate 
rough agreement between the computations, although more extensive calculations should be made, ft appears that 
the predicted pKp widths are always smaller (ranging from a factor of 2 to a factor of fO) than those of the 3 Pq model. 
Finally comparing 3 Pq and pKp to (the rather sparse) data indicates that the 3 Pq model tends to over-predict decay 
rates while the pKp model under-predicts them, ft is clear that more theoretical and experimental effort is required 
before this situation can be clarified. 

Finally, we compare the 3 Po model to all well measured rates in the PDG in Fig. [21 3 . As can be seen, the 
computation of strong widths is not as reliable as the computation of masses; a typical width error is 30% and can 
reach factors of 2 or even 3. 



3 The decay modes are as follows. [1] bi — * uitt, [2] n 2 — » /21", [3] Kq — ► Kit, [4] p — > mr, [5] <j> — + KK, [6] n 2 — * pm [7] tt 2 — * K* K + cc, 
[8] 7T2 -> wp, [9] 0(1680) -» K*K + cc, [10] K* -» Kit, [11] K*' -> Kit, [12] K*' ~^_pK, [13] K*' ~> K*tt, [14] D*+ -> D°tt+, [15] 
V>(3770) -> DD, [16] f 2 -» tttt, [17] f 2 -* KK, [18] a 2 -+ pn, [19] a 2 -+ jjtt, [20] a 2 -> KK, [21] f 2 -» KK, [22] D s2 -* DK + D*K + D s n, 
[23] K 2 -> Kn, [24] K 2 -+ K*tt, [25] K 2 -> pK, [26] K 2 ujK, [27] p 3 -» mr, [28] p 3 -> lott, [29] p 3 -+ KK, [30] K 3 -+ pK, [31] 
K 3 -+ K*tt, [32] K 3 -+ Kit. 




(23) 



TABLE IV: Open-flavor Strong Decays of Selected Charmonia. Widths assigned in this review are denoted 
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FIG. 2: 3 Pq Couplings Required to Reproduce Data|2f| . 



4- Critique 



Constituent Quark Models 

Most complaints about constituent quark models centre on one theme: quark models are not QCD. But models 
do not seek to be QCD, rather they attempt to capture the dominant physics relevant to the problem at hand; in 
this case, the structure and dynamics of hadrons. In this regard, most quark models suffer from serious deficiencies 
which will be become apparent as states higher in the spectrum are discovered. Namely: (i) Fock sector mixing 
(or coupled channel effects) will become more important as more continuum thresholds open up. Of course virtual 
meson loops also impact the low lying spectrum, but these effects are more easily absorbed into model parameters |27| . 
(ii) Although gluons may be integrated out to good accuracy low in the spectrum, they must manifest themselves 
as dynamical degrees of freedom somewhere (roughly 1300 MeV above the lowest S-wave states in a given flavour 
sector), (iii) The use of chirally non-invariant effective interactions (for example, scalar confinement and the 3 Pq 
model) cannot be accurate for systems with typical momentum scales well above the chiral symmetry breaking and 
quark mass scales |28j. 

Short Range Interactions 

Although one may hope that short range quark interactions are described by one gluon exchange, this is far from 
obvious in the intermediate regime relevant to hadrodynamics. In particular, spin-dependence in the quark interaction 
may arise due to one gluon exchange, coupled channel effects, relativistic dynamics, and instantons. The last three 
remain serious challenges for model builders. 

A good example of the difficulties which remain to be surmounted is the use of the scalar confinement interaction 
discussed in Section H B II Among its many problems are that it breaks chiral symmetry (and therefore is only useful as 
an effective interaction low in the spectrum where explicit chiral symmetry breaking is relevant) and that it cannot be 
used to describe baryons (quark-quark scalar interactions are of opposite sign to quark-antiquark scalar interactions, 
implying that baryons are dominated by an unphysical negative long range interaction). Indeed the entire analysis 
of Section TlB II was much too naive: QCD does not generate quark interactions of the structure J JKJ. In fact, 
spin-dependent interactions arise from nonperturbative gluodynamics involving virtual hybrid states. This can lead 
to complex behaviour - including permitting a vector confinement potential to have an effective scalar spin-dependent 
interaction. How this can be achieved is illustrated in Ref. • 

Lattice 

Although lattice gauge theory is 'in principle' equivalent to QCD, in practice it has far to go. A simple way 
to see this is by counting parameters - 'in principle' there should only be one parameter necessary for all lattice 
computations, a s , or equivalently, Aqcd- In practice there are many hidden parameters, including the lattice spacing, 
the lattice volume, the structure of the interpolating field (which typically must be carefully tuned to obtain a signal), 
Monte Carlo parameters (number of sweeps, Metropolis parameters, etc), autocorrelation parameters, chiral tuning 
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parameters, plateau parameters, etc. Of course the benefit of the lattice is that all of these dependences may be 
categorised as systematic errors and eventually reduced; something which cannot be said of models. 

However, more serious concerns remain. The computations reported in Table ITU employed the quenched approxima- 
tion, which unfortunately is not an approximation, but a truncation which renders QCD a sick theory. Furthermore, 
disconnected diagrams are rarely included in spectrum computations. Although we look forward to the day when these 
issues are overcome a warning is appropriate: unquenched lattice computations permit hadron decay, both virtual 
and real. Accommodating these in a rigorous way can only be achieved by computing observables such as scattering 
matrix elements. This represents a fearsome technical challenge. It also implies that full lattice computations will 
require all of the techniques of (and be subject to all of the ambiguities associated with) experimental partial wave 
analysis. Finally, it is worth stating that quenched lattice computations tend to succeed where constituent quark 
models succeed and fail where they fail. The convergence is not surprising: quark models fail when states have 
a complicated Fock structure and the lattice has difficulty with such states because they represent multiple scale 
unquenched systems. 

Decay Models 

The importance of continuum channels on the spectrum and hadronic properties has already been mentioned. Of 
course attempts to describe such physics relies on a thorough understanding of the transition operator which mixes 
Fock sectors. As discussed above, our understanding of these operators is essentially nonexistent. This leads to a 
variety of issues when attempting to 'unquench' the quark model. 

First, as opposed to quark model lore, continuum thresholds affect states throughout the spectrum. However, these 
effects may be largely renormalised away leaving residual effects near thresholds. Of course, the problem in hadronic 
physics is that continuum channels tend to get dense above threshold. Summing over these channels is nontrivial 4 - 
indeed the sum may not converge. 

Furthermore, one expects that when the continuum virtuality is much greater than Aqcd quark- hadron duality 
will be applicable and the sum over hadronic channels should evolve into perturbative quark loop corrections to the 
Wilson loop or quark model potential. Correctly incorporating this into constituent quark models requires marrying 
QCD renormalisation with effective models and is not a simple task. 

Finally, pion and multipion loops can be expected to dominate the virtual continuum component of hadronic states 
(where allowed) due to the light pion mass. This raises the issue of correctly incorporating chiral dynamics into 
unquenched quark models. The relationship of chiral symmetry breaking to the constituent quark model has been 
discussed in Ref. |3dj| and a variety of hadronic models which incorporate chiral symmetry breaking exist |3l| but 
much remains to be achieved. 



4 Note that the Cornell group only sum over six continuum channels, DD, DD* , D* D* , D a D 3 , D S D*, D*D*; the sensitivity of their 
results to this truncation should be tested. 
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II. THE X(3872) 

The A(3872) is the poster boy of the new heavy hadrons - it has been observed by four experiments in three 
decay and two production channels and continues to refuse to fit into our expectations for charmonium. Experimental 
properties of the A will be reviewed in the next subsection. Subsequent subsections will present model and effective 
field theory attempts at understanding the X. Finally, other heavy quark molecular states are considered. 

A. Experiment 

A summary of the discovery modes of the A (3872) is presented in Table IV! 



TABLE V: A(3872) Discovery Modes. 



mass 


width 


production/decay mode 


events 


significance experiment 


3872.0 ±0.6 ±0.5 


< 2.3 90% C.L. 


B ± -> A*X -> A^+tt" 


J/ij> 25.6 ± 6.8 


10(7 


Belle[32] 


3871.3 ±0.7 ±0.4 


resolution 


pp — > X — » TV + 1T~ J/lf) 


730 ± 90 


11.6a 


CDFII[33] 


M(J/ip) + 774.9 ± 3.1 ± 3.0 


resolution 


pp — > X — > TV + 1T~ J/ip 


522 ± 100 


5.2a 


D0f34] 


3873.4 ± 1.4 




B~ -> K~ X -> A-tt+tt" 


J/ip 25.4 ±8.7 


3.5a 


BaBar[35] 



The world average A mass is 

M(A) = 3871.9 ± 0.5 MeV (24) 

assuming independent systematic errors. 

The clarity of the A signal is exemplified in the irirJ/ip invariant mass distribution as seen by CDF-II and shown 
in Fig. 
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FIG. 3: The X in mrJ/ip^. 

As will be discussed below, the invariant mass ttt: distribution carries important information about A decays. Thus 
it is of interest that CDF-II have collected enough statistics to verify that the 7T7t invariant mass distribution in A 
decays is dominated by an intermediate p meson. This is illustrated in Fig . 0] 

BaBar[23 have set an upper limit on the r]J/%jj decay mode of the A|37l|: 

Br{B ± -> AA^) Br(X r]J/i()) < 7.7 ■ 10" 6 90%C.L. (25) 
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CDF II Preliminary, 360 pb"' 
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FIG. 4: 7T7r Invariant Mass Distribution 36 



The relevant data are shown in Fig. [3] 




JA|»r| Mass(GeV/cT) 
FIG. 5: Search for X -> nJlib^\. 

If the X is an isovector state charged partners should exist and one expects that the branching ratio for the charged 
X should be twice that of the neutral X . BaBar 381 have searched for these and find the limits 



Br(B° -> K + X~) Br(X~ -> ir~Tr°J/tp) < 5.4 • 10 



-6 



and 



Br(B- K°X~) Br(X- -► vr^vr J/i/>) < 22 • 10 



-6 



(26) 



(27) 



Spurred on by tetraquark models of the X (see Section Hi Bl . BaBar 39| have examined the mass difference of the 
Xs produced in charged and neutral B decays and determine this to be 



AM (X) = 2.7 ± 1.3 ± 0.2 MeV. 



(28) 



They have also determined that the ratio of neutral to charged production of the X is constrained to lie between 15% 
and 134%: 



Br(B° -> K Q TnrJ/ib) 
at 90% C.L. This result has recently been improved by BaBarj3^, who measure 5 



(29) 



Br(B° -> XK°) 
Br(B+ -> XK+) 



0.50 ±0.30 ±0.05. 



(30) 



5 This result has been updated from the previous 0.61 ± 0.36 ± 0.06 l40l . 
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Motivated by the predictions of Ref. Belle have searched for and found the X in the channel 7J/ , 0^^ thus 

establishing that the X has positive charge conjugation parity. The signal is shown in Fig. EJand has a significance 
of 4 a. 
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FIG. 6: Observation of X -> yj/ip^. 

Belle also have measured the product of branching fractions 



Br(B -> KX) Br{X -> jj/tp) = (1.8 ± 0.6 ± 0.1) • 10" 



and the ratio of widths 



(31) 



T(X -> 7J/» 

T(X -> TT+TT-J/V') 



= 0.14 ±0.05. 



(32) 



The same paper reports the discovery of the X decaying to ir + ir tt q J/iIj with a significance of 4 a (see Fig. [7J with 
a strength that is comparable to that of the ttttJ/i/j mode: 



Br(X -> tt+tt-tt J/i/') 



= 1.0 ±0.4 ±0.3. 



Br(X -> vr+vr- J/i/-) 

This observation establishes strong isospin violation effects, which will be discussed further in Section III Bl 
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FIG. 7: Observation of X -> 7r7r7r,7M|4^. 
The products of branching ratios have been measured by BelleQ and BaBar[39l IZF 



Br-(B + -> XK + )Br{X -> ttttJ/i})) = (1.3 ± 0.3) • 10" 



and 



(34) 



Br(B + -> XK + )Br(X -> rrrrJ/i/;) = (1.01 ± 0.25 ± 0.10) • 10" 



(35) 



15 



Finally, there are rumours that Belle observes the X decaying to D°D°ir at a rate ten times larger than that of 
mrJ/ip^^. The new decay mode is seen in 11.3 ± 3.6 events at 5.6 a with a product of branching fractions of 



Br(B -> KX) Br(X -► DDn) = (2.2 ± 0.7 ± 0.4) • 10" 



(36) 



Determining the quantum numbers of the X is a high priority experimental task. BES 46( have searched for the X 
in e + e~ collisions and find the bound 



r(e+e~) Br(X -» ttttJ/t/-) < 10 eV (37) 

at 90% C.L. Thus the X is unlikely to be a vector charmonium state. This conclusion is confirmed by Babar 41] who 
search for the X in initial state radiation and find an upper limit to the product of branching ratios: 

T(X e + e~)Br(X ->■ ■KitJ/iji) < 6.2 eV (38) 

at 90% C.L. 

Belle have measured angular correlations in the irnJ/ip system and use these to rule out possible J PC = ++ or 
^ values for the X. Furthermore the dipion mass distribution strongly disfavours 1 or 2 h assignments for the 
X^3|. Thus, it is li kely that the X is a 1 ++ charmonium state |44l|. 

CDF |47j and D0|34| have examined the production characteristics of the X. In particular, it is found that approx- 
imately 16% of the X sample at CDF comes from B decay; the remainder comes from direct production or deca y of 
short-lived particles. The mix of production sources is similar to that observed for the '0(25') charmonium state |47| 
(see Fig. [SJ. Thus it appears that there is no penalty for producing the X in pp collisions. 
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FIG. 8: Comparison of event yield fractions for X(3872) and ip(2S). (a) PtUNi) > 15 GeV, (b) \y\(J/4>) < 1, (c) cos (9^ < 0.4, 
(d) effective proper decay length < 0.01 cm, (e) isolation-1, (f) cos# M < 0.4|34l|. 

In summary, a narrow state at 3871.9 MeV is seen in B decays and pp collisions decaying to irirJ/ip, mnrJ/ip, 



PC 



r 



jJ/ip, and D D n which has production characteristics similar to the tl){2S). The quantum numbers J 
are strongly preferred by the data. Various attempts to interpret this data are presented in the following sections. 



B. Models 



Charmonium 



The first place to seek an explanation of the X(3872) is evidently in the charmonium spectrum jig. A glance at 
Table ITT1 reveals that the only viable candidates for the X are 2P or ID states, with the x' c o being the most likely 



16 



candidate based solely on mass. Table llVl indicates that 3 Pq widths are too large to accommodate a narrow 2P or 
ID charmonium (however, the pKp model yields smaller widths). Subsequent experimental effort (described above) 
has focussed attention on 1 ++ charmonia, of which the x' c i 1S the only possible candidate. But the predicted width of 
this state is 127 MeV in the 3 Pq model and 16 MeV in the pKp model - much too large. The apparent lack of viable 
charmonium candidates has led to much theoretical speculation. 

Three options are available when the quark model fails us: (i) attempt to change the quark model; (ii) seek 
supernumerary states; (iii) claim experimental or interpretational error. The first option is extremely difficult, as 
the average mass errors reported in Table IIIII illustrate. In general one must find an effect which is localised in the 
spectrum so as not to destroy the good agreement seen elsewhere. 

In the list of usual supernumerary suspects are hybrids, glueballs, diquark clusters, and molecular states. All of 
these options have been suggested for the X. The prominent ones are discussed in the following. 

2. Tetraquarks 

Maiani et al.^j^ have constructed a model of the X(3872) (and several other states as we shall see) which postulates 
that it is dominated by a diquark-diquark structure. The results of a relativistic tetraquark model are reported in 
Ref. 50]. Chief assumptions in the model are that quarks group in colour triplet scalar and vector clusters and that 
the interactions of these are dominated by a simple spin-spin interaction 51] . Using the X(3872) as input, Maiani et 
al. predict the spectrum of new diquark charmonia shown in Fig. [5] 
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FIG. 9: Spectrum of Tetraquark Charmonium States |4?S|. 

This spectrum assumes that the mass of the scalar and vector [cq] diquarks are 1933 MeV and constructs states as 
follows: 

(39) 
(40) 

l) + |Hv[cg] s ;J = l)) (41) 

l)-\[cq] v [cq} s ;J = l)) (42) 

(43) 
(44) 

Thus the model predicts that charged partners of the X(3872) should exist, namely X + = [cu][cd\ and X~ = [cd][cu]. 
Furthermore, the two neutral X states comprised of [cm][cu] and [cc?][cd] should mix with angle 6 and hence produce 
two neutral Xs with a mass difference of 

M(X h ) - M(X t ) = 2{md '2 u) = MeV * ( 8 ± 3 ) MeV ( 45 ) 

cos(26>) cos(2#) 
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(similar considerations apply to all the states in Fig. |5J. Furthermore, one of these states will populate B + decay 
while the other will populate B° decay. 

Relative strengths of the production of tetraquarks may be simply estimated. However, I disagree with the argu- 
ments of Ref. 49] where it is claimed that the amplitude for a strange quark to combine with the spectator quark is 
comparable or less than the amplitude for it to combine with a produced quark. Rather this amplitude is suppressed 
by a factor of 1/N C due to its colour structure. Thus one concludes that B° — > K ~X + , B° — > K°Xd , B + — » K + X u , 
and B + — > K°X + should all occur at the same rate. 

The model makes strong predictions which may be compared with experiment. For example, the ratio of neutral 
to charged production of the X is predicted to be unity whereas BaBar measure (Eq. |3*UI) a ratio of 50% with large 
errors. Unfortunately, many of the other predictions of the model are not as successful. The predicted mass difference 
of the neutral Xs (Eq. I45f) is 2 a removed from the measured mass difference (Eq. I28fl and no charged X states, or 
any of the other states of Fig. [5J have been sighted. Finally, the existence of diquark correlations in hadrons implies 
an extensive new exotic and cryptoexotic spectroscopy for which very little evidence exists 6 . 

While the ideas surrounding diquarks and tetraquarks are interesting and even compelling 52] , one must reserve 
judgement on their reality until strong evidence for them comes to light. 

3. Additional Models 

Li has suggested that the X is a ccg hybrid meson[53j which decays predominantly via X — > J/ipgg — > J/^pmr and 
has a significant J/ipa branching fraction. The existence of the opposite G-parity 3n decay mode is explained as due 
to isospin violation. Its large magnitude (see Eq. I33fl is left as a conundrum. 

Current lattice [BH and flux tube model[55| expectations for the mass of ccg hybrids are in the range 4200-4400 
MeV, thus there must exist significant systematic error in these models to permit a charmonium hybrid at 3872 MeV. 
In the absence of countervailing data, one must therefore regard hybrid models of the X as unlikely. 

Seth[5^] has proposed that the X is a vector glueball with a small admixture of vector cc. This proposal has since 
been eliminated by the observation of the uiJ/ip (Fig. 0) and jJ/tp (Fig. EJ decay modes. 

On the basis of a simple chromomagnetic interaction model, H0gaasen et aL[57j postulate that the X is a loose four 
quark agglomeration which is a mixture of cucu and eded flavour states. Thus another neutral X is expected along 
with two charged states cued and edeu. In addition, a broad isospin multiplet is predicted at 3742 MeV. 

Kalashnikova has constructed a simple coupled channel model of the cc spectrum which incorporates a central 
confining potential and 3 P model couplings to DD, DD*, D*D* , D S D S , D S D* S , and D*D* continual . Couplings 
are computed assuming SHO wavefunctions for the mesons and a simple renormalisation procedure is carried out. 
An interesting result of this computation is that the renormalised x'd S1 ^ s a ^ roughly 4000 MeV, while a virtual 
bound state appears just above DD* threshold. It is natural to identify this with the V(3872), and in fact, such 
an identification will be indistinguishable from a very weakly bound DD* molecule (to be discussed in the next 
section). Although this approach suffers from all of the problems mentioned in Section ll B 41 and in some cases can 
ruin agreement with naive quark model results (such as vector meson leptonic widths), this interesting observation 
warrants further investigation. 

Finally, Bugg 59] has reminded us that cusps occur in amplitudes at thresholds, and these can manifest themselves 
as bumps in cross sections slightly above threshold. The proximity of the X to D°D*° threshold implies that the cusp 
scenario should be treated seriously. In this regard, we note that typical enhancements due to cusps are O(Aquc) m 
width, in contrast with the very narrow X . Furthermore, if the X is simply due to the opening of a threshold, then 
the reaction B+ -> K°D + D*° should exhibit a peak similar to the X channel B+ ->■ K+D°D*°. That it does not 
indicates that significant final state interactions are affecting B + — > K + D° D*° . And of course these interactions may 
be strong enough to create a hadronic resonance. Additional observations on distinguishing bound states from cusps 
can be found in Ref. [7!| . 

4- Molecular Interpretation of the X 

In principle nothing in QCD prevents the formation of nuclear-like bound states of mesons and speculation on the 
existence of such states dates back thirty year s l60ll6H. T he proximity of the X to D"D"' threshold immediately led 
to speculation that the X is a DD* resonance [621 l63ll64l l65l log . Indeed D°D°* threshold is at 



The fiascos associated with pentaquarks and baryonia will only be mentioned in this footnote. 
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D°D°* = 3871.2 ± 1.0 MeV (46) 

while the world average X mass is (Eq. E3J) 3872.0 ± 1.8 MeV, 0.7 ± 1.2 MeV higher than threshold. The positive 
binding energy (with large error) is disconcerting for molecular interpretations, but, is no problem for tetraquark, 
cusp, or virtual state interpretations. In any event, the binding energy for a putative molecule is very small and this 
forms the basis for an effective field theory approach to the X, which will be discussed in the following section. 

In general, if the X is dominated by a loose D°D°* component then it does not have good isospin which should give 
rise to a distinct decay pattern. The most likely bound state is an S-wave which implies that the quantum numbers 
of the X are J PC = 1 ++ . Mixing with nearby Xci states should induce weak isospin violating decays in them as well. 
Finally, since the mesons are weakly bound they decay as if they were free, thus the X should decay via the modes 
X -> D°D°Tr and X -> D°D°j with a relative rate of 62:38. 

Making more detailed predictions forces one to employ a specific model because lattice gauge theory is not sufficiently 
well developed to permit accurate examination of this state (although a recent attempt has appeared see Section 
II B 41 for additional comments) and too many unknown couplings are present in effective field theory descriptions of 
the X. 

The first step in constructing a molecular (or any) model of the X is to determine the relevant degrees of freedom. 
The D°D*° channel has already been identified as important. This in turn implies that the binding energy is very 
small which implies that cc components in the X are small (this is quantified below). Similarly excluding hybrid 
components leaves meson- meson channels. Of these the four closest are listed in Table fVTI 



TABLE VI: Continuum Channels near 3872 MeV. 



channel 


threshold 


T(A) T(B) 




3871.2 


< 2.1 MeV 


D + D*~ 


3879.3 


96 ± 22 keV 


p°J/^ 


3867.9 


150.3 ± 1.6 MeV 91.0 ± 3.2 keV 


U}J/lp 


3879.5 


8.49 ± 0.08 MeV 



The next element in a microscopic description is dynamics. Two candidates come to mind: pion exchange and 
quark-level interactions. Indeed, if quark dynamics are limited to quark/colour exchange interactions then it is 
natural to admit both interactions since pion exchange is long range while quark exchange is a short range effect. 

The quark model employed here assumes nonrelativistic quark dynamics mediated by an instantaneous confining 
interaction and a short range spin-dependent interaction motivated by one gluon exchange. The colour structure is 
taken to be the quadratic form of perturbation theory. This is an important assumption for multiquark dynamics 
which has received support from recent lattice computations for both confinement |68| and multiquark interactions |69| . 
The final form of the interaction is thus taken to be 



— - — i"- i^ s '- s -{—)' 7' (47) 

where A is a colour Gell-Mann matrix, a s is the strong coupling constant, b is the string tension, rrii and rrij are the 
interacting quark or antiquark masses, and a is a range parameter in a regulated spin-spin hyperfine interaction. The 
parameters used were a s = 0.59, b = 0.162 GeV 2 , a = 0.9 GeV, and 0.335, 0.55, and 1.6 GeV for up, strange, and 
charm quark masses respectively. Relevant meson masses obtained from this model are p = 0.773 GeV, J/ip = 3.076 
GeV, D = 1.869 GeV, and D* = 2.018 GeV, in good agreement with experiment. 

Meson-meson interactions are obtained by computing the Born order scattering amplitude shown in Fig. 1101 for 
a given process [T^. Because of the colour factors in Eq. 07| this amplitude necessarily involves an exchange of 
quarks between the interacting mesons. Thus the leading order DD* interaction couples DD* with hidden charm 
states such as pJ/4> a- n d loJ/^j. This amplitude may be unitarised by extracting an effective potential and iterating 
it in a Schrodinger equation|7(j. The method has been successfully applied to a variety of processes such as KN 
scattering [jl} and J/tp reactions relevant to RHIC physics jT^. It has even proven surprisingly useful for relativistic 
(and chiral) reactions such as 7T7t scattering [7fjj. 

An effective potential is extracted by equating the scattering amplitude to that obtained for point-like mesons 
interacting via an arbitrary S-wave potential. The resulting potential is shown in Fig. 1111 The distinctive "mermaid 
potential" seen here is due to destructive interference between diagrams in the quark level amplitude. Thus details of 
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FIG. 10: Quark Exchange Diagrams Contributing to Meson-Meson Scattering. 




FIG. 11: Effective Potential for DD* — > ujj/tp. The dashed lines indicate different partitions of the four-quark Hamiltonian. 

the potential are sensitive to the assumed microscopic interaction, however, its general shape and strength are quite 
robust. 

Long range interactions are assumed to be dominated by one-pion-exchange. Ref. [(5(| chose to follow the method of 
Tornavist|6l| in constructing an effective pion- induced interaction. This is based on a microscopic quark-pion inter- 
action familiar from nuclear physics (in more modern language, the form of the interaction is dictated by spontaneous 
chiral symmetry breaking): 



L = - 



V2U 



Sx^x^^^d^^x). 



(48) 



Here f„ = 92 MeV is the pion decay constant, r is an SU(2) flavour generator, and g is a coupling to be determined. 
The effective potential is derived by projecting the quark level interactions onto hadronic states in the nonrelativistic 
limit. In the case of pseudoscalar- vector states one obtains [6lJ (see Appendix B for more details): 



-7V0 



1 
1 



C(r) 



V2 



-V2 



T(r) 



(49) 



where 



C(r) 



H 2 e'^ 
ml m-rrT ' 



T{r) = C(r) I 1 + — + ] 



(50) 
(51) 
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and 

3 9 
777 /7 

^"^l" 1 - 3 ^- (52) 

The matrix elements refer to S- and D-wave components of the pseudoscalar- vector state in analogy with the deuteron. 
The strength of the interaction has been fixed by comparing to the irNN coupling constant via the relationship 
/wAr/^ 71 " = 25/18 • m^g 2 /f%. This allows a prediction of the D* decay width which is in good agreement with 
experiment 61] . The parameter \x is typically the pion mass, however, one can incorporate recoil effects in the 
potential by setting p 2 = m 2 — (my — m p s) 2 ■ Finally, the coupling 7 is a spin-flavour matrix element which takes 
on the following values: 7 = 3 for I = 0, C = +; 7 = 1 for I = 1, C = — ; 7 = — 1 for I — 1, C = +; and 7 = — 3 
for I = 0, C = — . Thus the isoscalar positive charge parity channel is the most likely to form bound states and 
subsequent discussion focusses on it. 

The potential of Eq. |^]is often regulated to represent the effect of additional unknown short distance interactions. 
The regulator scale, A can be fixed by comparison with nuclear physics; for example NN interactions yield preferred 
values for A in the range 0.8 GeV to 1.5 GeV depending on model details. Alternatively, reproducing the deuteron 
binding energy requires A w 0.8 GeV. Tornqvist has employed an intermediate value of A = 1.2 GeV which is 
appropriate for D mesons and this is taken as the canonical cutoff in the following. 

Unfortunately, D mesons are sufficiently light that the DD* system (just barely) does not bind with canonical 
parameters. However, the combined pion and quark induced effective interactions are sufficient to cause binding. 
Thus it is of interest to study the properties of possible bound states as a function of their binding energy. This has 
been achieved by allowing the regulator scale to vary between 1.2 and 2.3 GeV. Binding is seen to occur for A larger 
than approximately 1.23 GeV. 

Wavefunction probabilities are defined by 



Z a = I \ip a \ 2 (53) 

where a is a channel index. These are shown as a function of binding energy in Fig. ^| It is clear that the D°D°* 
component dominates the wavefunction, especially near threshold. However, the D + D~* component rises rapidly 
in strength with isospin symmetry being recovered at surprisingly small binding energies (on the order of 30 MeV). 
Alternatively, the cuJ/ip component peaks at roughly 17% at Eb ~ 9 MeV. The contribution of the pJ /tp wavefunction 
remains small, peaking at less than 1% very close to threshold. 

100 
80 
60 

"a 

0- 40 
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FIG. 12: Wavefunction Components of the v|g(|. 




At the time of publication this model had two post-dictions: a charmonium resonance of mass 3872 MeV with a 
decay mode ofirirJ/ip- Predictions included the quantum numbers J PC = 1 ++ ; that the tttt invariant mass distribution 
should reflect its origins as a virtual p; and a 3nJ/ip decay mode of strength comparable to the 2ir mode with an 
invariant mass distribution dominated by the uj. All of these predictions were subsequently confirmed experimentally 
(see Figs 0] [7| and Eq. l3lfl) . Decays are discussed in Sect. Ill Dl Further data and discussion of the molecular 
interpretation of the X appears in Sect. Ill El 
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C. Effective Field Theory 

The notion that the X is dominated by a weakly bound D°D*° system suggests that an effective field theory 
approach to the X may be useful, in analogy to the successful description of the deuteron[73|. The utility of this 
approach is further enhanced by an approximate heavy quark symmetry and the anomalously small binding energy 
of E B = -0.5 ± 0.9 MeV, which should be compared to the natural scale of m„/(2fj, DD *) ps 20 MeV|73.l75j. 

Weakly bound systems may be described in terms of a single parameter, the scattering length a - a property known 
as 'low energy universality'. Simple results which apply to a weakly bound system include estimates of binding energy 



and the weak binding wavefunction 



ip(r) 



-r I a 



(55) 



or 



(56) 



The scattering lengths can be quite large (compare to the deuteron, cld = 4.3 fm). For example, Eb = 0.5 MeV 
corresponds to a = 6.4 fm, while Eb = 0.1 MeV implies that a = 14.4 fm 7 . Universality ideas were first exploited by 
Voloshin[74[ who estimated the effects of interference on the processes X — > D°D°j and X — > D q D°tt°. 

Braaten and Kusunoki have followed this initial work with a series of papers examining production and decay of 
the X. This began with the observation that the probabilities for non-D°D*° components of the X, X production, 
and X decays all scale as 1/a and hence are strongly suppressed in the weak binding limit. 

The analysis was performed with the aid of an effective field theory which couples D, D* and \ mesons j?^: 



n 



m D o (D^D + D^D) + m D *o (D f • D + D f • D) 

+ (m D a + m D *o + ^o)^ 1 ■ X 



1 



(D^V 2 D + D^V 2 D) 



1 



2m£)0 2mj3*o 
-\{m D a +m D ,o)- 1 x t • V 2 x 

A (DI) + L»D) f • (DD + GD) 

+.9o [x f • (DD + DD) + (DD +DD) t - x ] 



(D f • V 2 D + D f • V 2 D) 



(57) 



The bare x binding energy is represented by vo- The four body interaction of strength Ao represents short range 
interactions between D and D* mesons while the term proportional to go couples the x to its constituent DD* 
channel via a pointlikc interaction. 

Low energy universality also permits analysis of 'coalescence' models of X production in T and B decays (see 
Fig. H3J with the results^ 



r(T(45) -> X + hadrons) - E B \ogE B 



and 



A(B -> XK) 



1/2 



2ckp B ■ (*x 



(58) 



(59) 



7 A subtlety associated with these scales will be discussed below. 
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FIG. 13: X Production from T(4S) Decays via Coalescence (left)[7(| and from B Decays (right) |7^|. 

for X production in B decays[72}- Here A is an ultraviolet cutoff, and c is a coupling denned by A(B — > DD*K) = 
cpB ■ e. This expression was analysed further in Ref. [T^ where the authors set the cutoff to be A = w m w which 
is a 'cross over momentum' scale defined by the transition from a phase space dominated invariant mass distribution 
to a resonance dominated invariant mass distribution. The improved estimate is then 



Braaten and Kusunoki also obtained 



Br(B° -> XK°) , . 

Br { B + ^XK+) = -° m - < 61 > 

(The upper end of this estimate agrees with the lower end of the molecular estimate to appear in Eq. 1711 ) 

Finally, Braaten and Kusunoki have used factorisation between universal l ong range behaviour and unknown short 
range behaviour to predict the binding energy dependence of decay rates [79ll80| . Using the results of Ref. (6|| to fix 
the short range results, yields predictions such as 



T(X -> W/0 = 222 keV (^J^) ^ • (62) 

This method describes decays such as X — > irirJ/ijj via a virtual pJ/ip channel and associates that channel with the 
scattering length a. It should be noted, however, that the presence of multiple channels generates multiple scattering 
lengths via the coupled channel weak binding relations 



E B (a) = T (63) 

where a denotes an X channel and Es(a) = m to t{a) — mx- Thus, although the probability of channel a scales as 
l/o, the dynamics associated with that channel scale with a a = 1/ y/2fi a EB(ot). 

These scales can be seen clearly in the X molecule model of the previous section. Numerically integrating the 
coupled channel Schrodinger equation and fitting the subsequent channel wavefunctions to the weak binding form 
(an example is shown in Fig. 114ft yields estimates of the channel scattering lengths. For example, setting A = 1.13 
GeV gives a binding energy of 0.9 MeV and hence dpofiro — 23.8 GeV" 1 while the fit to the numerically obtained 
wavefunction yields a D ojj,o = 24 GeV -1 . Similarly, for the pJ/ip wavefunction one obtains a fit value of a pi p = 9.5 
GeV -1 while using Eb{p^) — 9.2 MeV gives an expected value of a p ^ = 9.3 GeV -1 . Similarly a D - D *+ = 7.6 GeV -1 . 
A selection of X properties as a function of A are given in Table EU Thus many scales are present when multiple 
channels contribute to the molecular structure of the X. As the binding energy decreases these additional channels 
disappear as l/a D o£,*o and the scale associated with each channel saturates at a a = [2/i Q (m tot (a) — m tot (_D°l)* ))] -1 / 2 . 
At larger binding energies the interplay of these scales can complicate the description of X dynamics and care must 
be exercised. 

Notice that the weak binding wavefunction of Fig. 1141 is quite accurate for distances greater than 2 fm. Similarly, 
the Fourier transformed weak binding wavefunction is accurate for momenta less than 200 MeV. Above this scale the 
wavefunction is quite small, hence the weak binding wavefunction may serve as a useful approximation to the full 
wavefunction in favourable circumstances. 
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TABLE VII: X Properties vs. Cutoff 



A (GeV) 


E R (MeV) 




^ I 






1.14 


1.0 


89% 


4.0% 


5.0% 


0.7% 


1.17 


1.1 


86% 


5.5% 


6.4% 


0.8% 


1.20 


1.4 


82% 


7.5% 


8.1% 


0.8% 


1.23 


1.7 


78% 


9.2% 


9.3% 


0.8% 


1.30 


2.7 


69% 


14% 


12.7% 


0.7% 


1.40 


4.8 


60% 


20% 


15% 


0.5% 


1.50 


8.0 


55% 


24% 


16% 


0.3 % 



0.3 
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o 
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0.1 
0.05 


10 20 30 40 50 60 70 80 90 100 
R (GeV 1 ) 

FIG. 14: The full DqDq wavefunction and the weak binding wavefunction. 




Finally, Al Fiky et al. [8l| have improved the effective field theory of Braaten and Kusunoki (Eq. I57|l by incorporating 
heavy quark and chiral symmetries. Thus the theory is written in terms of a superfield doublet combining the 
pscudoscalar and vector mesons: 



H. 



(Q) - 



1+ j6 



H 



<«>« = T °£T(«)t 7 °. 



(64) 



These fields have the usual transformation properties under heavy-quark spin symmetry and SU(2)y flavor symmetry 
and describe heavy mesons with a definite velocity v. One then has C = £2 + £4 with 



£2 = - «Tr 


H (Q) v ■ VHW 


- 1 Tr 








2nip 





ig^ — (Q) 



(65) 



and 



-U = ^Tr 



7m 



Tr 



^Tr 

4 



Tr 



H^m (Q) r l5 



(66) 



The covariant derivative is given by T>£ b — Sabd^ — (1/2) [yd^ 1 ^ + £d^£>) , and g is the D*Dir coupling. The pion is 

realised nonlinearly as £ = exp(i7r • r/y/2f„). The third term in Eq. (|65|l is needed to account for the D — D* mass 
difference M(D*) - M(D) = -2\ 2 /M(D)_ 

Evaluating the traces yields for the DD* sector 



C4,dd> = - C 1 D«)t jD (Q) £ ,*(Q)t £ ,*(Q)M _ Cl£) *(Q)t jD *(Q)M jD (Q)t jD (Q) 



(67) 
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This lagrangian differs from that of Eq. 1571 where the interaction is described in terms of a single parameter Ao- 
This is related to the parameters above by Ao = —C\ = C2. The additional parameter in the lagrangian of Al Fiky et 
al. may be ascribed to distinguishing effective pseudoscalar and scalar exchange interactions between heavy mesons. 
Nevertheless, the binding energy of heavy pseudoscalar-vector states depends on a single linear combination of C\ 
and C2. The mass of the X may be used to fix this linear combination and a prediction for a BB* bound state at 
10604 MeV is obtained 81]. Predictions for possible DD, BB, BD* , and DD* bound states are not possible until 
another linear combination of the couplings can be determined. Further discussion of other molecules is contained in 

sect, mm 

The mass of the BB* molecule was obtained with the aid of the heavy quark mass scaling relation Cj ~ 1/toq which 
implies that Ci{rrib) = m c /mbCi(m c ). Furthermore, the pion exchange portion of the coupling constant (obtained 
upon integrating out the pions) is much smaller than the contact portion|82|. Both of these observations are not in 
agreement with the microscopic model quark+pion exchange model: pion exchange dominates the effective interaction 
and the pointlike interaction is controlled by scales such as the string tension, the light quark mass, and the pion 
mass. 

The power and limitations of the effective field theory approach are easily seen. When enough data exists to 
determine coupling constants reliable predictions are possible in a broad range of problems. These predictions are 
under-pinned by QCD symmetries or universality. An example is the scaling of decay rates with binding energy 
shown in Eq. 1621 However, the method fails when data is sufficiently sparse that parameters cannot be determined 
or multiple scales affect the dynamics. 

D. Microscopic Decay Models 

1. Strong Decays 

The difficulties associated with predicting strong decays of mesons have already been discussed in Section H B 41 
These difficulties are worse in molecular decays where uncertainties are compounded by poorly constrained structure 
and dynamics. There is, however, no utility in being too easily dissuaded! In this spirit three- and two-body decay 
processes driven by one gluon exchange are displayed in Figs. 1151 and ITol respectively. The latter replace constituent 
mesons with kinematically allowed decay product mesons, either by rearrangement or by flavour changing processes. 
The former proceed via quark pair creation such as drives 3 Pq or Cornell model meson decays. In the case of 
the X molecular model of the previous section, dissociation gives rise to decays such as X — + pJ/ip — > mrJ/ip or 
X -> DD* -► DD-y. 




FIG. 15: X Strong Decay via Dissociation (left) and mr Production (right). 




FIG. 16: X Strong Decay via Rearrangement (left) and Annihilation (right). 



In terms of the molecular X model, the leading source of decays is dissociation. In the weak binding limit the 
decay rates are simple to estimate: to good approximation dissociation will proceed via the free space decay of the 
constituent mesons. Thus, for example, 



T(X — > TTlTj/ljj) = Z plf ,r(p — > 7T7T) 



(68) 
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and 

T(X -» D°D°tt°) = 2Z 0Q T(D*° -> D°n°). (69) 

The factor of two in the latter expression accounts for the S-wave combination -±=(D°D*° + D°D*°) in the bound 
state wavefunction. A more extensive list of predictions is given in Tabic IVIIII 



TABLE VIII: Some Decay Modes of the X(3872) (keV). 



Be (MeV) 




D°D° 1 


D+D-n 


(D + D°-k-+c.c)/V2 £T 


"L>"7 


7T+7T J/tp 


TV + TV yJ/lp TV + TT 7r°J/%p 


7T°7J/V> 


0.7 


67 


38 


5.1 


4.7 


0.2 


1290 


12.9 720 


70 


1.0 


66 


36 


6.4 


5.8 


0.3 


1215 


12.1 820 


80 



Strong decays in the case of strong binding are considerably more complicated. The formalism which describes this 
is discussed in Appendix A. 

2. Radiative Decays 

Radiative decays of molecules provide an interesting contrast to the impulse approximation processes which drive 
mesonic radiative decays. In the case of the X radiative decay can occur via vector meson dominance (in the uJ/ip 
and pJ/tp components) or annihilation (DD* components), as illustrated in Fig. El The vector meson dominance 
diagram will contribute to X — > ^yj/ip. It is clear that the amplitude must be proportional to the light vector meson 
wavefunction at the origin and to the \ wavefunction for the channel in question. The specific result is 

Tvmd = ^f^\Mr = 0)| 2 + 3Z H 2 Ml)) 2 ■ (™) 

Decays to jip' or 7^2 can only proceed via annihilation and hence are very small. 




FIG. 17: X Radiative Decay via Vector Meson Dominance (left) and Annihilation (right) [4^ . 

A variety of predictions are presented in Table 17x1 Here, the fourth column summarises the results of Barnes and 
Godfrev|48| which are computed in the impulse, nonrelativistic, zero recoil, and dipole approximations. The results 
for X — + jJ/ip are particularly sensitive to model details - in fact this rate is identically zero for SHO wavefunctions. 
The sensitivity is examined in columns five and six which present the results of two additional computations. The 
first, labelled [A], employs the same approximations of Barnes and Godfrey but uses meson wavefunctions computed 
with a simple, but accurate, Coulomb+linear+smeared hyperfine potential. It is apparent that the 7 J/ip rate is very 
sensitive to wavefunction details. Furthermore, one may legitimately question the use of the zero recoil and dipole 
approximations for the jJ/tp mode since the photon momentum is so large in this case. The sixth column (model 
[B]) dispenses with these approximations, and one finds a relatively large effect for "fj/ip. 

Tabic IIXI makes it clear that computations of the 7 J/ip radiative transition of the x'ci are very sensitive to model 
details 8 . The result of Barnes and Godfrey is similar to that computed here for the molecular x state but is much 
smaller than models A and B. Furthermore, the rates for 7^ and 7^2 are very small for a molecular X (at the 



There is an additional error induced by arbitrarily changing the quark model \' mass to 3872 MeV. This is made clear through the 
observation that the dipole formula for the width scales as g 3 whereas the momentum space formula scales as q. 
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order of eV) and quite small for a charmonium X. Perhaps the most robust diagnostic is the yip' decay mode. For 
a molecular \ this can only proceed via the annihilation diagram and hence is very small. Clearly a measurement of 
the yJ/tp and ytp' decay modes of the AT (3872) will provide compelling clues to its internal structure. 



TABLE IX: El Decays of the X(3872)|42 



mode m f (MeV) q (MeV) 



r[cc] (keV) T[cc] (keV) T[cc] (keV) 
[B&G] [A] [B] 



r[Xci] (keV) 



yj/tp 
7^'(2 3 Si) 
7^"(l 3 Di) 



3097 
3686 
3770 
3838 



697 
182 
101 

34 




8 
0.03 





E. New Data and Further Analysis 



The apparent success of the molecular interpretation of the X has been challenged by several recent experimental 
results. These are reviewed here. 

• The ratio of production of the X in neutral and charged B decays (Eq. I30[) lies in the range 0.25 - 0.97 with 64% 
C.L. and hence is rather large compared to expectations based on the molecular picture. This ratio may be estimated 
with the diagrams shown in Fig. ^| The colour leading diagram (left panel) permits the processes B + — > K + D°D* Q 
which will contribute to K+X, B+ -> K°D a D* + which tests X+ production and cusp effects, B° -> K°D-D*+ 
which accesses the X through the suppressed charged component, and B° — > K + D~ D*° which tests X~ production 
and cusp effects. The colour suppressed diagram (right panel) permits the following decays: B + — ► K + D° D*° , 
B+ -> K+D+D*-, B° -> K°D°D*°, and B° -> K°D + D*~ . Since the diagrams involve identical mesons (at the 
strong level) and only differ by topology one may predict the ratio 



Br(B° -> K°X) _ 
Br(B+ -> K+X) 



,1/2,2 



I4Z, 



1/2 
00 



Z_ 



1/2 1 



0.06- 0.29. 



(71) 



The lower estimate assumes that there is no relative phase between the topologies (it is difficult to imagine why there 
should be one) and that the X is completely dominated by the neutral DD* component. The higher estimate assumes 
that the neutral component is 90% of the X while the charged is 10% 9 . Given the large errors in the experimental 
result of Eq. [201 it is too early to rule out the molecular interpretation of this state. An accurate determination of 
this ratio will measure the relative importance of the suppressed D + D*~ component of the X. 





FIG. 18: B — > KX. Colour enhanced amplitude (left). Colour suppressed amplitude (right). 
Eas. 1531 [551 and I5B1 imply that 



Br(X -> D a D°7T a ) 
R = r i t^t = 22 ± 13. 

Br(X -> 7T+7T- J/iP) 

This should be contrasted with an earlier measurement from Belle which determined 



(72) 



9 Assuming a negative relative phase raises this estimate to 0.25 - 0.51. 
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Br{B+ -> XK+)Br(X -► D°D ir°) < 6 • 10~ 5 (73) 

at 90% C.L. Eqs. and then imply that i? < 6. A reliable determination of this ratio will be crucial to 
understanding the X and further experimental effort is encouraged. 

In the molecular picture R is driven by decays D*° — > D°ir° and p — * irir and hence is given (in the weak binding 
limit) by 



2Z oa T(D*°) 

R - z*m • (74) 

With the pion and quark dynamics assumed in Section III B 41 one obtains Zqo rs 0.8 and Z p ^, w 0.01 at a binding 
energy of 0.7 MeV. Thus R « 0.08, clearly in conflict with experiment. 

If the X is indeed a weakly bound molecule, Eq. EJ] must be accurate and hence the only way to reconcile this 
prediction with the data is if the pJ/tp component has been severely overestimated. The required probability is 
Z p ^k. 0.007%. 

• The large rate for D Q D°ir° implies that Br(X -> ttttJ/^) < 0.05. With Eqs. 01 EU this implies that 



Br(B -> XX) > 2 • 10" 4 (75) 

which is comparable to the branching ratio for Xd- Br(B -±_XdK) = 4.0(1.0) • 10 -4 ^}- However, both of these 
results are very close to bounds established by BaBar |39l l40t l84j : 



Br{B ± -> K ± X) < 3.2 • 10~ 4 (90%C.L.) (76) 

and 

Br{X -> ttttJ/^) > 4.2%. (77) 

These results, and the conflict with the earlier Belle result (Eq. I73|) . can be most easily interpreteted by assuming 
that the rate to D°D°ir has been overestimated in Eq. E21 

If a large D°D°tt° decay mode is confirmed, the simplest interpretation is that the X contains substantial cc - 
in keeping with the production characteristics discussed at the end of Sect. Ill Al More specifically the X must mix 
strongly with the x' cl . This may be estimated with the aid of the 3 Po model and Fig. The net result is a mixing 
matrix element 



a x = V2Z^ 2 J d 3 kip x {k)A(-k) (78) 

where A is the 3 P amplitude for x'd — * D°D*°. Here tfjx(k) = {n y /a)^ 1 (k 2 + a -2 ) -1 is the weak binding X 
wavefunction. 




FIG. 19: X - Xd Mixing. 



Evaluating the integral yields surprisingly large mixing matrix elements, which are summarised in Tabled As can 
be seen, the mixing amplitude is diminished as the X scattering length gets larger. Surprisingly, mixing with the 
excited x'ci i s n °t strongly suppressed with respect to the ground state Xd- Most surprisingly, the mixing matrix 
element can be quite large, so large that the relatively close Xci(3930) has a very (nonperturbatively) large mixing 
amplitude with the X. This unexpected result explains the large KX branching ratio in B decays and the X 
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TABLE X: X 


- Xci 


Mixing. 




state E B (MeV) a (fm) 


Zoo 


a x (MeV) 


prob 


Xci 0.1 14.4 


93% 


94 


5% 


0.5 6.4 


83% 


120 


10% 


Xci 0.1 14.4 


93% 


60 


100% 


0.5 6.4 


83% 


80 


> 100% 



production characteristics in pp collisions. The complete phenomenological implications of this result have not been 
explored [85j. 

• Belle have found the jJ/ip decay mode of the X and report a large branching of 14% of the rnrJ/ip mode (see 
Eq. 1321) . In the model of Sect. Ill Dl one expects this to be 0.4% at a binding energy of 0.7 MeV - much too small to 
explain the data. Note, however, that the predicted rate for x'ci ~* lJ li> is 10 - 140 keV (Table llX|l ; the upper end 
of this prediction implies a branching ratio of 14% if the TnrJ/ip mode has a width of roughly 1 MeV. But this implies 
that the X has a very large x'ci component or that the predicted irirJ/ip rate is too large. In view of the previous 
discussion, both of these possibilities appear likely. 

F. Summary 

It is possible that all of the new experimental data may be explained if the coupling of the X to the x'ci 1S correctly 
incorporated in the coupled channel formalism. This would increase the jJ/ip branching fraction, permit the large 
D a D a ir a decay mode, and allow the charmonium-like production characteristics. A detailed phenomenology of this 
scenario remains to be constructed. 

Note that cusp, tetraquark, renormalised x'cii virtual state, and molecular interpretations all assume that the X 
is dominantly (or heavily modified by) a ccuu Fock space component. The chief difference among the models is the 
spatial, colour, and spin configuration of the state, which is of course, driven by different assumed dynamics. Thus 
distinguishing these models is a qualitative issue and choosing between them is question of efficiency of description 
and consistency of adjunct predictions. 

Finally, we comment on two misapprehensions in X dynamics that have appeared in the literature. First, the 
SttJ/i/j decay mode proceeds via the ujj/ip component of the X. It is sometimes thought that the decay occurs via 
the natural line width of the to. Since the width of u> is only 8 MeV and the wJ/ip system is 7.5 MeV above the X 
mass, the probability of the X(u>J/ip) fluctuating into ZirJ/ip is very small. Thus the large relative branching fraction 
for the 3tt mode cannot be correct. This argument is incorrect: the ojJ/ip component of the X is off-shell due to the 
dynamics leading to binding and the u need not fluctuate to a "lower mass" to permit decay. 

Second, Suzuki[8(| has argued that the one-pion-exchange potential as described in Sect. Ill B 41 is not valid because 
the process D* — > Dn can occur on-shell. This argument is incorrect. The four point function that describes DD* 
scattering does indeed exhibit a pole (as opposed to the NN interaction) due to set of measure zero in the phase space 
available to the DDir system. The pole induces a small imaginary part to the scattering amplitude, which describes 
the decay X — > DDit. Of course, this decay may be separately estimated in perturbation theory as in Sect. Ill Dl 
The real part of the scattering amplitude remains unchanged and its analysis may proceed as for NN scattering. The 
procedure is described in more detail in Appendix B. 

G. Other Molecules 

Other heavy molecular states can be generated from the one pion exchange mechanism discussed above. Effective 
potentials in the various channels are discussed in Appendix B. Isovector channels never bind; some properties of 
isoscalar channels which bind are listed in Table |X"T1 These results have been computed with A = 1.23 GeV and 
generally agree with T6rnqvist[6lJ except that I find only one D*D* bound state with J PC = 0++. 

Note that binding energies in the BB channels are quite large. In particular, the l ++ (BB*) has a binding energy of 
43 MeV, in contrast to the binding energy of 0.2 MeV predicted by Al Fiky et al.^,^ (see Sect. ITTT^ . The latter 
is predicated upon an effective field theory description of HH* interactions and the scaling relation C\ ~ 1/toq. This 
is quite different from one pion exchange and quark exchange interactions which do not scale with tuq. Thus there 
is a conflict between these approaches (where in principle, none should exist) and further effort is required to resolve 
this discrepancy. Certainly, one expects that the BB* system should be more strongly bound than DD* since the 
kinetic energy is greatly reduced - in agreement with the pion exchange model predictions. 



29 



TABLE XI: Isoscalar Heavy Quark Molecules 



tit fifp 






mass CMpVI Rb 


D*D* 


0++ 




4019 1.0 


BB* 


o-+ 


3 Po 


10543 61 


BB* 


1++ 


3 S U 3 L»i 


10561 43 


B*B* 


0++ 


1 5o, 5 -Do 


10579 71 


B*B* 


o-+ 


3 Po 


10588 62 


B*B* 


1+- 


3 Si, 3 £>i r 


10606 44 


B*B* 


2++ 


1 D2, 5 <S I 2, 5 £>2, 5 G2 


10600 50 



The molecular Xco (the J PC = ++ D*D* state) is interesting because it may be accessible in B factories and at 
Fermilab. This state can decay to wJ/tp, r)r) c , and r)'r] c . A rough computation of the partial widths yields 



T(Xc0 -» ujJ/ip) w 200 keV 



(79) 



and 



r(x c o - Wc) « 600 keV. 



(80) 



Finally, unlike the case of the X, no nearby hidden flavour meson-meson channels exist for these states. Thus their 
short range structure is expected to be substantially simpler than that of the X(3872). 
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III. OTHER NEW CHARMONIA 



The h c has been observed by CLEO[87| in the reaction ip(2S) 



(77) (iVc) with a mass of 3524.4 ± 0.6 ± 0.4 



MeV at a significance greater than 5 a . The observation has been made in exclusive measurements where the r\ c has 
been reconstructed, and in inclusive measurements (see Fig I20[) . The mass may be compared to the quark model 
predictions listed in Table ITT1 which range from 3516 to 3526 MeV (or 3474 MeV for quenched lattice computations). 
Furthermore, the measured product of branching fractions 

Br(^(2S) -> n°h c )Br(h c -> 7 ry c ) = (4.0 ± 0.8 ± 0.7) • 10~ 4 (81) 

is in close agreement with a theoretical estimate of 3.8 • 10~ 4 of Ref. |88J|. (This reference may also be consulted for an 
extensive discussion of spin splittings in quark models.) Thus it appears that the h c is conforming to expectations. 



400f 




3.50 3.52 



3.54 3.50 3.52 
M(jt°recoil)(GeV) 



3.54 



FIG. 20: Inclusive Production of the h c in 4>(2S) -* n°X. Monte Carlo (left); Data (right) |g2 



Nevertheless, as discussed in Section II B II the h c is an important diagnostic for the spin-dependent structure of 
the quark-quark interaction. The fact that it lies within 2 MeV of the spin-averaged mass of the \c.J multiplet: 
M c .o.g. — ^(xo + 5xi + 9x2) = 3525.36 implies that the long range spin-spin interaction V4 is very small, in agreement 
with the effective scalar confinement picture. 



As discussed above, hyperfine splittings are a useful diagnostic of spin-dependent confinement interactions. An 
old result from Crystal Ball 89] caused some consternation since it implied a rather large spin splitting in the radial 
ip(2S) — r] c (2S) system. The situation has recently changed with new measurements from Belle, CLEO, and BaBar. 

Belle 90] have observed an rj' c candidate in 39 ± 11 events in the KsKtt system of the decay B — ► KKsK~n + . The 
measured mass of the resonance is 3654 ± 6 ± 8 MeV and the width was determined to be T = 15^11 MeV ( r < 55 
MeV at 90% C.L.). 

Belle also measured the ratio of branching fractions 



Br(B -> Krj' c )Br(ri' c -> K s K--k + ) 
Br(B -> Kr) c )Br(r) c -> K s K~n+) 



= 0.38 ±±0.12 ±0.05. 



(82) 



Assuming equal branching ratios for rj c and rj' c to KsKtt (this assumption is discussed in Ref. [9ll ]'l gives Br{B 
Kr)' c ) I Br(B — > Krj c ) — 0.38. This may be compared to the factorisation expression: 



Br(B -> Krf c ) _ P/ (rf c ) 



f 2 . 



(tub + m K ) 2 — m 



Br{B->Kr, c ) p f (n e ) tf c |e(w^ c ))| 2 [(m B + m K f 



(83) 



v . v '= (m| + m 2 K - mt 1 )/(2m B m K ). 



where w 

A quark model computation gives [9 
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f Vc = 424 MeV 



and 



f v f e = 243 MeV. 

And, using a simple model for the Isgur-Wise form factor £(u>) = (2/(1 + w)) 2 gives 



(84) 
(85) 



Br(B -> KrQ 
Br(B -> #r? c ) 



0.38 



(86) 



in agreement with the data and the assumed equality of branching ratios. 

More recently, the CLEO collaboration 93] have observed the v( c in 77 — > 77^ — > KsK^tt^ (Fig. right panel) and 
find a Breit-Wigner mass and width of 3642.9 ± 3.1 ± 1.5 MeV and T = 6.3 ± 12.4 ± 4.0 MeV respectively. Assuming 
equal KsKir branching fractions as above, they also obtain 



r(r?c -> 77) = 1.3 ±0.6 keV. 



(87) 



Finally, BaBar[94| have observed an r]' c candidate in 112 ± 24 events in 77 — > rj' c — > KsKir (Fig. [2] left panel). 
They fit a mass of 3630.8 ± 3.4 ± 1.0 MeV and obtain a width of T = 17.0 ± 8.3 ± 2.5 MeV. The angular distribution 
of the decay products and of the total transverse momentum both suggest an assignment of J PC = v or J > 2, 
with the latter being unlikely for such a light state. 
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FIG. 21: Fit r) c {2S) by BaBar|ii|(left) and CLEO|a3l(right). 

As mentioned above, the rj' c is of some interest because its splitting with the is driven by the hyperfine interaction 
and hence probes this interaction in a new region. Specifically, the ground state vector-pseudoscalar splitting is 



M(J/i>) - M(r? c ) = 117 MeV 
whereas the excited splitting is now measured to be 



M(ijj') - M{rf c ) = 48 MeV. 



(89) 



Theoretical expectations for the former range from 108 to 123 MeV in simple (or 'relativised') quark models (see 
Table [nj and thus are within expectations. Alternatively, the latter is predicted to be 67 MeV in the quark model 
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of Eichtcn, Lane, and Quigg 95] . However, the authors note that including unquenching effects due to open charm 
meson loops lowers this splitting to 46 MeV which is taken as evidence in favour of their 'unquenched' quark model. 
However, the simple quark models mentioned above^^ find splittings of 42 - 53 MeV, indicating that it is too early 
to make definitive conclusions about loop effects. It is worth noting, furthermore, that attempts to unquench the 
quark model are fraught with technical difficulty (see Section ll B 41) and a great deal of effort is required before we 
can be confident in the results of any model. 

C. X(3940) 

The Jf(3940) is seen by Belle [96|j recoiling against J/if) in e + e~ collisions. The state has a Breit-Wigner mass of 
3943 ± 6 ± 6 MeV and a width of less than 52 MeV at 90% C.L.0. The X is seen to decay to DD* and not to uij/ip 
or DD. 

It is natural to attempt a 2P cc assignment for this state since the expected mass of the 2 3 Pj multiplet is 3850 
- 3980 MeV (see Table IHJ) and the expected widths are 20 - 130 MeV (see Table HVJ). Indeed, if the DD* mode is 
dominant it suggests that the V(3940) is the x' c i- There is, however, a problem with this assignment. Fig. 1221 shows 
the recoil mass spectrum up to 4.5 GeV. The bumps visible in the figure are a 3.3 a fluctuation at 2500 MeV 10 , the 
i] c at 2980 MeV, a structure usually associated with the x c o(3415) (but probably due to Xco, Xa, an d Xci)-, the Vc 
at 3638, and the V(3940). Thus there is no clear evidence for the Xci, and one may suspect that the x'd should 
also not be a strong signal in this reaction. This has led to speculation that the V(3940) is the radially excited 
77" 98, 99]. Unfortunately this interpretation also has its problems as the expected mass of the 77" is 4040 - 4060 MeV, 
approximately 100 MeV too high. An angular analysis of the decay products of the X will help resolve this issue. 




FIG. 22: V(3940) Signal in e+e~ -> J/ip+ charm. Ref. @. 

An interesting addendum to the issues surrounding the V(3940) concerns attempts to compute the production 
of the various charmonium states. Table IXLfl displays the predicted cross sections from two NRQCD computations 
and experimental results from Belle and BaBar. It appears that N RQCD is under-predicting r\ c J ftp production cross 
sections by an order of magnitude. However Bondar and Chernyak 103] have argued that NRQCD is not applicable 
to processes such as e + e~ — > J/ipr] c because the charm quark is not sufficiently heavy. Using general methods for 
hard exclusive processes and a model J/ tp l ight cone wavefunction gives the result labelled as BC in Table IXIII The 
last line refers to a light cone computation |l04| which also produces much larger cross sections than NRQCD. The 
authors also not that NLO corrections can be large. 

TABLE XII: Cross Sections (fb) for e+e" -» J/ipH at ^fs = 10.6 GeV. 



Could this be the tensor glueball? 
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Vc 


D .» u,,,. Ton! 


1 / .0 ±1 Z.O ± Z.L 


lU.o rt Z.O it l.c 


ID. 4 ± O.I ± O.U 


BelleflOO] 


25.6 ±2.8 ±3.4 


6.4 ± 1.7 ± 1.0 


16.5 ±3.0 ±2.4 


BL[101] 


2.31 ± 1.09 


2.28 ± 1.03 


0.96 ±0.45 


LHCfl02] 


5.5 


6.9 


3.7 


BCfl03l 
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BLU104] 


26.7 
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D. F(3940) 

Belle claim the discovery |l 05) of a second resonance at 3940 MeV, this time in the w,J/ip subsystem of the process 
B — ► KiririrJ/t/j. The S-wave Breit Wigner mass and width are (see Fig. I23[) determined to be 3943 ± 11 ± 13 MeV 
and 87 ± 22 ± 26 MeV respectively. A total of 58 ± 11 events are claimed for the Y with a significance of 8 a. The 
state has not been seen in the decay modes Y — > DD or DD*. Because it does not share the production or decay 
channels of the X (3940) it is possible that these states are distinct. 

As for the X(3940), the mass and width of the Y suggest a radially excited P-wave charmonium. However, the 
wJ/ij) decay mode is peculiar. In more detail, Belle measure 

Br(B -> KY)Br(Y uoJ/ij)) = (7.1 ± 1.3 ± 3.1) • 10~ 5 . (90) 

One expects that the branching ratio of a radially excited \cj is less than that of the ground state \ c j\ Br(B — > 
Kx' cJ ) < Br(B — > Kxcj) = 4(1) • 10~ 4 [I3|. Thus Br(Y — > ajj/ip) > 17%, which is unusual for a canonical cc state 
above open charm threshold. 

Thus the Y is something of an enigma, driving the observation of the Belle collaboration that it is perhaps a 
hybrid. The unusual decay mode and branching fracti on ar e consistent with the claims that a hybrid should have 
strongly suppressed DD, DD* and D*D* decay modes [l06| . However, a mass of less than 4000 MeV is in conflict 
with quenched lattice computations of the low lying hybrid spectrum|54j. Thus more data are required before the 
exotic nature of the Y can be claimed with any certainty. 




3880 4080 4280 3880 4080 4280 
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FIG. 23: Threshold Enhancement in the ujj/ip System (left) and the Breit Wigner Fit (right) |l05j. 



E. Z(3930) 

The Z(3930) was observed by the Belle collaboration [lO^ in 77 -> DD with a mass of 3931 ± 4 ± 2 MeV and a 
width of 20 ± 8 ± 3 MeV at a claimed significance of 5.5 a (see Fig. l24|) . The DD helicity distribution was determined 
to be consistent with J = 2. Due to its mass and decay mode, the Z may be considered a x(2P) candidate along 
with the X(3940) and F(3940). In this case identification as the x' C 2 i s natural. 
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The predicted mass of the x' C 2 is 3970 - 3980 MeV (Table© and the predicted width is 60 MeV (Table IE). 
However, setting the mass to the measured 3931 MeV restricts phase space sufficiently that the predicted strong 
width drops to 35 MeV, reasonably close to the measurement. The predicted branching fraction to DD is 70%. 
Furthermore, the largest radiative transition is x' C 2 ~ * V>'7 with a rate of 180 ± 30 keV[l3T|. 

If this assignment is correct it implies a 40 MeV error in the quark model predictions, somewhat large compared 
to the average errors shown in Table Hill and this is perhaps an indication that the expected breakdown in the model 
is occurring. Furthermore, if the V(3940) is the x'ci then these excited P waves exhibit an inverted spin splitting of 
-10 MeV rather than the expected +45 MeV. The perturbative Xc2 — Xci splitting is given by (16a s /(5r 3 ) — b/r) /m 2 , 
so an inverted splitting is possible, but indicates a serious lack of understanding of the relevant dynamics. If this 
scenario is confirmed then the x' c q is expected at M(x' c i) ~ l/2(8a s /r 3 — b/r) /m 2 . This is difficult to estimate in view 
of the uncertainty in the dynamics. A reasonable guess may be a x'co mass of 3945 MeV. Alternatively, if the mass 
scale has only been misunderstood, then a x'co is expected at approximately 3810 MeV - it is possible that a glimpse 
of this state is already visible in Fig. 1241 

In view of the uncertainty in the masses of the x' c j it is worth examining the partial widths of the x'co and x' C 2 as a 
function of the resonance mass. The results are shown in Fig. 1251 and indicate that the x' c o can t> e very narrow near 
3920 MeV. In this case it may not be a surprise that this state is not strongly visible in Fig. 
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FIG. 25: DD and DD* Partial Widths vs. 



XcJ 



Mass. 



The production rate in 77 has been determined to be 



T 77 Br(Z -► DD) = 0.23 ± 0.06 ± 0.04 keV. (91) 

The 3 Po model results pr edict a DD branching ratio of 70% for the x' c2 and hence 1X77) = 0.33 keV. This rate has 
been estimated by Barnes |l08j to be 0.64 keV - given the uncertainties, this is reasonably close to the data. 

Searching for DD* decay modes of the Z will be an important diagnostic for this state. Measuring helicity angle 



35 



distributions for the A, Y, and Z will also be important for determining quantum numbers. Finally, increasing 
statistics for all these states will be useful. 



F. F(4260) 

The Y"(4260) is an intriguing recent addition to the charmonium spectrum. It was discovered by BaBar [l09j as 
an enhancement in the nnJ/ip subsystem of the reaction e + e~ — * JisrJ / ipmr (see Fig. The measured mass is 

4259 ± 8 ± 4 MeV and the width is 88 ± 23 ± 5 MeV; with a total of 125 ± 23 events in the signal. The rnr mass 
distribution peaks near 1 GeV (see Fig. I77|) . suggesting an isobar Y — > J/^/n (980) — > J/^tttt decay (although there is 
substantial strength at lower invariant mass suggesting a a contribution). The reaction rate has also been determined 
to be 



r(Y -> e+e - ) Br{Y -> J/^tttt) = 5.5 ± 1.0±g f eV. (92) 

The Y is not seen in e + e~ — > hadrons, in particular in the final states DD or D S D S , indicating that the branching 
fraction of the Y to J/ifiirir is large compared to that for the ^(3 770). 

BaBar have also searched for the Y in B decays and report [llCj "an excess of events" around 4.3 GeV at a 
significance of 3 a. The measured branching ratio is 



Br(B- -> K~Y)Br(Y -> J/^tttt) = (2.0 ± 0.7 ± 0.2) • 10" 5 . (93) 

If the J/ipmr mode dominates the width of the Y this implies a B decay branching ratio of approximately 10~ 5 , 
two orders of magnitude smaller that the branching fractions for KJ/ip or Kip'. Similarly, T(Y — > e + e~) ~ 5 eV, 
substantially smaller than typical vector leptonic widths. 




m(jtirJA|/) (GeV/c 2 ) 
FIG. 26: The Y(4260) Observed in e+e" -> -vr.g g .7M7r7r[Toflj. 

Evidently the Y is a vector with cc flavour content. Of course the low lying charmonium vectors are well known 
(see Table ITTf) and the first mesonic charmonium vector available is the 3 3 Z?i state. However quark model estimates 
of its mass place it at 4460 MeV, much too heav y for the Y. Of course this statement relies on the robustness of the 
quark model. Thus, for example, Llanes-Estrada flTH has argued that the Y is the ^/ , (4S I ) based on the spectrum of 
a relativistic model. This displaces the ^(4415) which moves to the ip(5S) position. 

A variety of o ther explanations fo r th e Y have been proposed. Among these are the suggestion that the Y is a 
A C A C barvonium |ll2j . Maiani et q/. |ll3| claim the Y is the first orbital excitation of a tetraquark [cs]g[cs]s state 
which decays predominantly to D S D S . A nonet of vector Y states is predicted with masses Y(I = 0) = 3910 MeV 
and Y(I = 1/2) = 4100 MeV. Alternatively, Liu, Zeng, and Li suggest that the Y is a XciP molecule bound by a 
exchange |l!4j . Thus two isotriplet partner states are predicted and possible isotriplet XcoP and Xc2P states may exist 
as well. 

Of the charmonium states which we know must exist, the most natural explanation is as a cc 
A hybrid assignment agrees with the large J/ipnir width 11 , small leptonic coupling, and small B decay branching 



Although it remains to be established that a large J/tpTnr mode is inconsistent with canonical charmonium. 
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FIG. 27: 7T7T Invariant Mass Distribution [123]. The line corresponds to an S-wave phase space model. 



fraction. As stated in the Introduction, the lightest charmonium hybrid is expected at 4400 MeV - somewhat high, 
but perhaps acceptable given our lack of experience in this sector. The observed decay to a closed charm final state 
is simply expl ained : all hybrid decay models forbid decays to identical mesons, and many forbid decays to pairs of 
S-wave mesons [lO^ j. Thus the allowed channels for a vector hybrid are DD* (which is suppressed), D^ p (a suppressed 
D-weive mode), DiD, and D[D. The latter modes are above threshold and thus are very small. Ref. |l!7j reminds us 
of these facts in the context of a constituent glue model 12 . Thus the dec ays o f a putative vector hybrid at 4260 MeV 
will be dominated by closed charm channels. The UKQCD collaboration |ll9j has found a large coupling of bb hybrids 
to closed flavour modes such as Xb®- Such a large coupling will be required to accommodate the relatively large total 
width of the Y. Understanding the mechanisms behind this coupling will then be a major task for phcnomcnologists. 

If the Y is indeed a vector cc hybrid, it should also decay to J/ijjcr, J/ipKK, Xco^, h c rj, and h c rf and these modes 
should be sought. Furthermore, one may e xpect nearby pseudoscalar and exotic J PC = 1 ^ hybrids which will 
display their own unique decay patterns |116| . 

The proximity of the Y to DD\ threshold (at 4285 MeV) suggests that the resonance could actually be a cusp 
at DD\ threshold combined with rescattering into the J/ipfo ~> J/V' 7r7r channel. Confirming the phase motion 
of the Y is thus an important (although very difficult) task. Searching for a DDttit decay mode with the Dtttt 
subsystem dominated by a virtual D\ would also be an important diagnostic. Furthermore, if rescattering drives the 
Y enhancement then one may expect similar enhancements at DD S \ (and similar) thresholds. These may possibly be 
seen in the J /ipKK channel. 

In line with the previous discussion, it is tempting to speculate on molecular interpretations for this state. In 
particular the DD\ channel is an S-wave threshold and is thus a likely one for forming a molecule. However, if the 
system does bind, it does so with a novel mechanism since pion exchange does not lead to a diagonal interaction in 
this channel (unlike the case of the AT (3872)). It is possible that off-diagonal interactions may provide the required 
novelty. For example, DD\ — > D*D* — » DD\ via iterated pion exchange; however, it seems unlikely that such a 
mechanism can lead to binding. Countering this, is a rece nt lat tice computation that claims to see a four-quark state 
at 4238 ± 31 MeV and predicts a cess at 4450 ± 100 MeV|l2(i|. 

Lastly, the reader may be d isturbed that the Y appears in ISR but is not seen e + e _ annihilation. This issue has 
been considered in Ref. |109| | where the Y cross section in e + e~ was estimated to be 50 pb. This corresponds to 
0.34% of the signal at -y/i = 4.25 GeV, which is not experimentally discernible. Thus the ISR production mechanism 
serves as a useful filter for some vector channels, such as ttkJ /4>- 



Detailed predictions of constituent glue models of hybrids should be regarded with care since they do not agree with level orderings of 
adiabatic hybrid surfaces as computed on the lattice unless specific gluodynamics are assumed lllSl . 
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IV. NEW OPEN CHARM STATES 



The new r evolu tion in charmonia was contemporaneous with BaBar's discovery of an enigmatic open charm state, 
the £>cf2317) |l2l| . We shall proceed with a review of theoretical expectations for heavy-light mesons before moving 
on to the properties and interpretations of the new open charm states. 



A. Heavy Quark Symmetry and the Quark Model 

1. Notation 

Generic Qq states will be denoted as D in the following. Properties of the six lowest lying Ds are listed in Table 
IXIIII Note that charge conjugation is no longer a symmetry in the heavy- light sector; thus the 3 Pi and 1 P\ states 
mix to form the physical D\ and D[. Quark model assignments for the remaining states are as in the quarkonium 
sector. 

In the uiQ — > oo limit the heavy quark acts as a static colour source and it is useful to describe the system in the 
jj coupling scheme. Thus the light quark degrees of freedom are written as j q — £ q + S q where i q (S q ) is the angular 

momentum (spin) of the light quark. The total spin-parity is then given by J p — | (8 j q - 

The heavy-light system can be very different from quarkonium. For example, quarkonia contain a hierarchy of 
scales, toq, (p) — a s mQ, E = oqrriQ, etc, whereas a heavy-light system has two scales, tuq and Aqcd- This can have 
profound implications, for example, hyperfine splitting for heavy quarkonium must scale as 8Mh vp ~ toq whereas it 
scales as 8Mh yp ~ ^■QCD/ m Q f° r heavy-light systems. 



TABLE XIII: Low-Lying Qq States 



property 


D D* 


D Di D[ D 2 


J p 

(2 S+ l) Lj 

£q 

Jq 

decay 


(T 1" 



1 - i - 

2 2 


0+ 1+ 1+ 2+ 

3 P 3 Pl - 1 Pl 3 Pl ~ X Pl 3 P2 
11 11 
1+ 1+ 3+ 3 + 

{DP) S (D 2 *P) S (D*P) D (D 2 P) D 



2. Spin Splitting and Multiplet Structure 

The different bases impose different multiplet structures on the D spectrum. In the absence of spin splitting, the 
spectroscopic LS basis has an S-wave (D,D*) doublet and a P-wave quadruplet. Alternatively, the heavy quark basis 
contains three doublets labelled by = \ ,\ , and | . 

The spin splittings induced in the LS basis may be determined by evaluating matrix elements of the interactions of 
Eq. |21 In the following only terms of order 1 /mt or 1/ m q m,Q are retained. The perturbative hyperfine splittings is 

SM(D) = (94) 
4 SmQm q 



5M(D*) = +\-^- (95) 
4 imqm q 



The P-waves are split by tensor and spin-orbit interactions according to 



SM{D ) = — — 



)>K, 



SM(D 2 ) = I {Vs ° 



2 



(96) 



(97) 
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Here Vso is given by 

V50 = f + ^ + ^ 5 (m 9 ,m Q ). (99) 

The sign of the off-diagonal term in the P-wave splitting is convention dependent. Here I have chosen to work with 
qQ mesons and to combine angular momentum and spin in the order L x S. 
Tensor splittings are necessarily of order l/mQm q and are given by 

6M(D ) = -4 (Fs) (100) 

llmQlTlq 

5M{D 2 ) = -\-P^L- (101) 

5 lZlJlQlTlq 

SMPPi) = +2 ^ (102) 
12toqto 9 

Diagonalising the spin-orbit splittings in the heavy quark limit yields 

M(D 2 ) = M(D[) = Mo + (103) 

and 

M(A)) = M(Di) = Mo - (104) 
in agreement with the jj multiplets. The physical D\ and 13^ states arc defined in terms of a mixing angle: 

l-Di) =+cos( ( /»)rP 1 )+sin(0)| 3 P 1 ) 

\D[) =-sin( ( /»)| 1 P 1 ) +cos(^)| 3 P 1 ). (105) 

In the heavy-quark limit a particular mixing angle follows from the quark mass dependence of the spin-orbit and 
tensor ter ms, w hich is 4>hq — —54.7° (35.3°) if the expectation of the heavy-quark spin-orbit interaction is positive 
(negative) |122| . Since the former implies that the 2 + state is greater in mass than the + state, and this agrees with 
experiment, we employ <f> = —54.7° in the following. This implies 

\Di)hq = +^Pi) -\/f |3Pl> 

\D[)hq = +y| 1 P 1 ) + -±=| 3 P 1 ). (106) 



3. Decays 



The P-wave D mesons are sufficiently light that the only isospin conserving decays are to DP or D* P where P is 
a pseudoscalar meson (jr for D decays and K for D s decays). Conservation of angular momentum implies that the 
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Dq can only couple to DP in an S-wave and one may thus expect the Dq to be broad. Similarly the D-i only couples 
to DP in a P-wave and is narrow. 

In general a D\ or D[ may decay to D* P in S- and D-wave. However, the conservation of heavy quark spin implies 
that the D\ only decays in an S-wave while the D[ decays in a D-wave. Thus the j q = | doublet is broad and the 

j q = I doublet is narrow. Fig. 1281 illustrates this selection rule |l22tll23j |. We consider an initial D^ decaying into a 
D* and a light pseudoscalar meson. Conservation of angular momentum and heavy quark spin implies that 



j q (DP)=j q (D*)+L 



D*P- 



(107) 



Setting j q (D*) = ^ and j q (D\) = \ then implies that the relative angular momentum in the final state is Lq'p = 0. 
Alternatively for the D[ one has j q (D[) — | and hence Ld*p = 2. 




FIG. 28: P-wave Decay in the Heavy Quark Limit. 

It is interesting to exam ine this behaviour in the 3 Pq model. One finds the following relation between 3 Pi and 1 P\ 
decay amplitudes [26l Il22] 



A{ l Pi -> VP) S 
AC Pi -> VP) D 



= -V=2 A{3Pl VP)i 
= V2A( 3 Pi^VP) D 



Eq. I105l then implies 



A{D 1 

A(D 1 - 
A(D[ 

A(D[ 



■VP) S 
VP) D 

■VP) S 
VP) D 



sin 6 1 A 



COSlj 

(V2 cos 4> + sin 4>)Ad 

/sin^ \ 
\^—j=-+cos cpJAs 

(— v2sin^!> + cos cj>) Ad 



(108) 

(109) 

(110) 
(111) 

(112) 



Substituting the heavy quark mixing angle cos <j> — 1/ v3 yields 

A{D 1 - VP) S = - 
A{D 1 -> 

.4^ -> ff) s 



2 






-4s 



(113) 

(114) 
(115) 
(116) 



in agreement with the heavy quark selection rule. The agreement is likely quite general because relations 1 1 081 hold in 
the Cornell model as well. This may be true because the nonrelativistic structure of the 3 Po and Cornell transition 
operators are very similar: 



H3P0 



(117) 
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while 

^Cornell ~ I • V(q) b\ +q a ■ qdl k . (118) 

k 

Additional selection rules may be derived|2^|. For example the following relations hold in the 3 Pq model 

A(2 3 S 1 ^ 1 P 1 P) S = -^=A(2 3 S 1 ^ 3 P 1 P) S (119) 

A(2 3 S 1 -> 1 PiP) D = V2A{2 3 S 1 ^ 3 P 1 P) D (120) 
A( 3 D 1 ^ 1 P 1 P) S = V2A( 3 D 1 ^ 3 P 1 P) S (121) 

Thus the Ditt (D[tt) mode will be large (small) in D(2 3 Si) decays and the Diir (D^ir) mode will be small (large) in 
L»( 3 L>i) decays. 

4- Spin Splitting with Continuum Mixing 

Although the conclusions reached above seem to be in accord with data (see Sect. IIV C(l . the logic used - namely 
first perform perturbation theory in the spin orbit interaction, then use the perturbed states to compute decay widths 
- is not correct. In general the coupling to continuum channels also perturbs the energy and structure of the -Dis; 
hence spin orbit and loop effects should be considered at the same time. In this case the perturbative mass splittings 
are given by 

SM(Do) = -=- + / ^—ttt (122) 



m 2 J (2tt) 3 Mo - E DP {k) 



m[D2) = + J w« (123) 

where uij(k) — (Dj\V\DP; fc) and V is, for example, the 3 P transition operator. 

The Di and D[ eigenstates are determined by diagonalising a Hamiltonian describing the interactions of the Pi, 
Pi, (DP)s and {DP)d channels. Integrating out the continuum channels and examining the perturbative mass shifts 
yields an effective 2x2 Hamiltonian to be solved for the physical Di states: 



SH 



1 (Vso) ( 1 -V2\ f d 3 k 1 / (2|^(fc)| 2 + |p 2 (fc)| 2 ) (-V2p (fc)| 2 + ^p 2 (fc)| 2 ) 

2 ml [-V2 ) + J (2n) 3 M -E D * P {(-V2\Q (k)\ 2 + ^ 2 (k)\2) (p (fc)| 2 + p 2 (fc)| 2 ) 
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where tDj,(fe) = ( 1 Pi\V\D*P; fc)^. The continuum contribution to 5H may be written as 



(124) 



STT - f d3k Po(fc)| 2 + P 2 (fc)| 2 (\ 0\ , f d 3 k p (fc)| 2 -lp 2 (fc)| 2 ( i _V2\ n9 ,, 
6M cont -J (2?r) 3 Mo _ EDtp ij+J^f M -E D * P [-V2 J [1Zi)) 

showing that the eigenvalues and mixing angle remain unchanged from the naive mixing case. 
The full perturbative mass shifts are 

m z J (2iry Mq- E D * P 



and 



*M(n>\ l <yso) * f d 3 k p 2 (fc)| 2 

SM{Dl) = 2^nf + 2j WMo-^-P (12?) 
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Finally, in the heavy quark limit Ep,p = Pd* p and ujo — — V3 Cjq. Thus the (Poj D{) doublet is recovered. However, 
0J2 = — -y/3/5 o>2 and the heavy doublet (D[, D 2 ) is not recovered. This is not a great surprise - while it was established 
that the heavy doublet was narrow, there is no reason to believe that the widths of the D[ and the D 2 are identical. 

Furthermore, it should be recalled that all possible virtual continua contribute to the mass splitting. These continua 
do not necessarily obey the amplitude relations given in Eq. 11081 and the naive heavy quark mixing angle and width 
relations need not hold. Also, the full computation should not be made in perturbation theory because some of the 
widths, and hence mass shifts, can be very large. In this case the concept of the mixing angle is not valid because the 
continuum cannot be represented as an effective potential but must be included in the Fock space expansion of the 
state 

|£>i) = Z l 1 /2 \ 1 P 1 ) + Z l 2 /2 \ 3 P 1 ) + J d 3 kZ{k)^ 2 \D*P:k). (128) 
Alternatively, the Hilbert space becomes infinite dimensional rather than simply two dimensional. 



5. D-waves 



Mixing in the 3 D 2 — 1 D 2 system may be addressed in the heavy quark limit of the constituent quark model as 
above. The mixing angle may be defined as 

\D* 2 ) = +cos(<j )D )\ 1 D 2 ) +sm( ( f >D )\ 3 D 2 ) 

\D* 2 ') = -sm{<t>D)\ 1 D 2 )+cos(4> D )\ 3 D 2 ) (129) 
Proceeding as with P-waves one finds a 3 D 2 — 1 D 2 mixing matrix of 

_ (Vso) ( 1 s/l 

Diagonalising the spin-orbit interaction yields the results 



M( 3 D 2 ^D 2 )=M a -^l ^ VD) (130) 



Af( 3 ^) = M(D* 2 ) = M - \{H q so ) D (131) 
M( 3 D 3 ) = M{D* 2 ')=M Q + {H q so ) D , 
and a 3 D 2 — X D 2 mixing angle of 4>d = —50.76°. Thus, the heavy quark P-wave states are 

2,lrvx ,/3 |3; 



\D* 2 )hq = \/ 7 /D 2 )-^-\ 3 D 2 ) 
v'\ .. - . / 3 |ln_\ 1 . / 2 |3; 



\D* 2 )hq = ^\ L D 2 ) + ^-\ 3 D 2 ). (132) 

As with P-waves, the 3 Pq strong decay model makes specific predictions for D-wave heavy quark decay amplitudes 
which may be useful in interpreting the spectroscopy. Some of these are: 

ACD 2 ^VP) p = -^A( 3 D 2 ^VP) P (133) 

A{ 1 D 2 ^VP) F = +^A( 3 D 2 ^VP) F (134) 

A{ l D 2 -» 1 P 1 P) = (135) 

The second of these is an example of the 3 Pq selection rule forbidding such transitions among qq spin singlets. 

We note that, in analogy with the P-waves, the amplitude ratios above imply that the D 2 decays strongly in P 
wave while the D 2 decays only in F wave and is thus narrower than the D 2 . As with P-waves, this conclusion agrees 
with spin conservation in the heavy quark limit. Unfortunately, the ability to distinguish the states is weakened by 
the many other decay modes which exist for these states. Nevertheless, the D 2 and D 2 have total widths which 
depend strongly on (pp> and measurement of any (or several) of the larger decay modes will provide (over) constrained 
tests of the model and measurements of the mixing angle. 
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B. Chiral Symmetry Breaking and Restoration in Heavy-Light Mesons 

1. Effective Field Theory 



Ten years prior to the discovery of the Z? s (2317), Nowak et aL |124| and Bardeen and Hill 125] investigated the 
consequences of an effective lagrangian description of heavy-light mesons, called the chiral doublet model. Nowak 
et al. assumed heavy quark symmetry and broken chiral symmetry and thus arrived at the lagrangian of Eq. 1651 
Alternatively, Bardeen and Hill assumed that chiral symmetry was manifest and thus wrote an analogous lagrangian 
with the linear realisation of chiral symmetry. The properties of heavy-light mesons have also b een e xamined in the 
context of a Nambu-Jona-Lasino model with chiral and heavy quark symmetries by Ebert et aL |l26| . 

The combination of heavy quark and chiral symmetries imply that the bare low-lying heavy-light states are degen- 
erate multiplets, H and G (or H') where 

ff=((T,l-) (136) 

and 

G=(0+1+). (137) 

Note that this axial is the j% = ^ + state. Heavy quark symmetry imposes the doublet structure shown, while chiral 
symmetry implies that they are degenerate. This structure is modified in two ways: the heavy quark symmetry is 
explicitly broken by the finite heavy quark mass and chiral symmetry is spontaneously broken by nonperturbative 
vacuum gluodynamics (explicit chiral symmetry breaking is neglected). The subsequent mass splittings raise the G 
multiplet above the H by roughly the chiral symmetry breaking scale and induce identical hyperfine splittings in the 
G and H multiplets. Thus one has the prediction 

M{Dx) - M(D*) = M(D ) - M(D) « A QCD . (138) 



Both sets of authors have subsequently applied the formalism to the newly discovered D s states 127]. As we shall 
see, Eq. 11381 is quite accurate: 

M(D sl ) - M(D*) = 347.2 ± 1.5 MeV, M(D s0 ) - M(D) = 349.1 ± 1.0 MeV. (139) 

Additional predictions are |l27| 

M(D ) = 2215 MeV (140) 

and 

M{D X ) = 2357 MeV. (141) 

Experimentally these are 2352 ± 50 MeV and 2421.8 ± 1.3 MeV respectively 0. 
In the B system the predictions are |l27| 

M{B ) = 5627 ± 35 MeV (142) 

and 

M{Bx) = 5674 ± 35 MeV. (143) 

The PDG reports several narrow and broad resonances of undetermined quantum numbers at 5698 MeV, in rough 
agreement with these predictions. 

The chiral nature of the lagrangian also permits computations of strong transitions involving Goldstone bosons. 
Indeed, parity doubling implies that the transition strength governing (0 + , 1 + ) — > (0~, l~)ir is related to the scalar- 
pseudoscalar mass difference by a Goldberger-Treiman relationship 
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ft,/* = M(0+) - Af (0-) (144) 

where f„ is the pion decay constant. Unfortunately, these decays are kinematically forbidden for the D s system and 
they can only proceed via isospin breaki ng D s — ► Drj —> Dir° transitions. Similarly couplings to multiple pions may 
be computed; for example Bardeen et aZ. |l27| predict 

T(D sl -> D s nn) = 4.2 keV. (145) 

It is worth stressing that the relationship shown in Eq. 11381 is very difficult to obtain in the constituent quark 
model. In particular, the Do and D\ are P-wave excitations with small hyperfine splittings. Thus the two sides of the 
equation are not related in any simple way. This issue will be examined in detail in the following section. 



2. An Explicit Chiral Doublet Model 



The speculation that chiral symmetry is relevant to heavy- light mesons is a vital part of the precedi ng ar guments. 
While general arguments have been made that this should be true high in the meson or baryon spectra |l28|, it is not 
clear that chiral symmetry is relevant to low lying states. Thus it is useful to examine the issue of chiral symmetry 
restoration in the meson spectrum with the aid of a simple relativistic Nambu-Jona-Lasinio quark model. The model 
is defined by the Hamiltonian 



H = J d 3 x^(-ia-d + /3m)^ + ^ J d 3 xd 3 y p a {x)V(x - y)p a (y) (146) 

where (3 is the Dirac matrix 70, a = ft^f, and p a is the quark charge density current. 

Chiral symmetry breaking occurs for a wide range of interactions, V and is embodied in the gap eauation |l29j| : 



C F 

Z^(A)ps p = [m(A) + Sm(A)\ c p + — 



d3k vi v\ 



(2tt) s 



■SfcCp p ■ k CfcSj 



(147) 



The functions and are defined in terms of the Bogoliubov angle 4>{k) as Sk = sin^(fc) and Ck = cos </>(£;). The 
constant and the UV cutoff A are required for renormalisation. However, this will not be discussed further since 
it is irrelevant to the following. 

The quark gap equation is to be solved for the unknown Bogoliubov angle, which then specifies the quark vacuum 
and the quark field mode expansion via spinors of the form 



Comparing the quark spinor to the canonical spinor (nonrelativistic normalisation is used) permits a simple interpre- 
tation of the Bogoliubov angle through the relationship p(k) = ktarnp(k) where p, may be interpreted as a dynamical 
momentum-dependent quark mass. Similarly p,(0) may be interpreted as a constituent quark mass. 

Chiral and heavy-quark symmetry can be studied by projecting the Hamiltonian on pseudoscalar, scalar, vector, 
and axial vector Qq channels. For example, the pseudoscalar and scalar Tamm-Dancoff Salpeter equations are 



and 



Mo-4>(k) = (E k + E k )<p(k) + l -J -0^ 

1 f d 3 p 



ap ■ k + b V(p — k) <fi(p) 



M 0+ 4(k) = (E k + E k )4>(k) + 



2 J (2tt) 3 



bp ■ k + a 



V(p-k) <t>(p). 



(149) 



(150) 



The factors in the interaction are given by 
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a = V 1 - -'W 1 _ -^V 1 + s 9^ /1 + s k + V 1 + s W l + s k\/l- v/l - s fc ( 151 ) 

and 

6 = v/1 + S,V1 + Sky/1 + s q Vl + s k + y/1 - s q y/l - s kx /l - s q Vl - s k . (152) 



The function s k is related to the dynamical quark mass via s k = A*(^)/ \J^ 1 + M 2 ^)- 

The structure of these equations simplifies substantially in the nonrelativistic limit, s k ,s k — * 1. In this case the 
pseudoscalar interaction reduces to 2 J V(j> while the scalar interaction reduces to 2 J p ■ kVtfi. Thus the pseudoscalar 
interaction is S-wave, while the scalar interaction is P-wave and one expects the scalar mass to be larger than the 
pseudoscalar. Of course these conclusions are in agreement with constituent quark models. 

Alternatively, these equations reduce to 

M ±<jy(k) = (E k + E k )ct>(k) + ( 44 \p-k+l] V(p - k) 4>{p) (153) 



(2tt) 3 

in the heavy quark (s k — > 1) and chiral (s k — > 0) limits. A similar results applies to the axial and vector equations. 
Thus the expected multiplet structure is recovered. 

The effects of spontaneous chiral symmetry breaking may be examined by setting s k = fi(k)/E(k). In this case the 
(Do, D\) doublet is split from the (D, D*) doublet. Hyperfine interactions arising from a finite heavy quark mass then 
split the doublets. Bardeen and Hill assume that this splitting is identical in both doublets, leading to Eq. 11381 This 
is consistent with the analysis here if the chiral splitting is small, or equivalently, (s k ) 1. However, experience with 
light mesons indicates that (s k ) « 1 (or (k) « fJ,(0)). If this is also true in heavy-light mesons (and we shall argue 
that it is) then the chirally split Dq and D are dominantly P- and S-wave respectively, and hence have very different 
hyperfine splittings and the conclusions of |l24l are incorrect. In general one expects (k) ~ /i(0) ~ A-qcd 

so that (sfe) = 0(1). Alternatively, the function fi is universal (it is determined by the structure of the vacuum) and 
hence the ratio 

<*>"*•- '"W 8 - (2)V»„i. (154) 



(k)\ p V fi(p) 

In summary, the scenario of Nowak et al. and Bardeen et al. relies on effective chiral symmetry restoration in 
low-lying heavy-light mesons but there is no reason to believe that this actually occurs. Of course this is a theoretical 
belief which flies in the face of the data (see the subsequent discussion); however, one must allow the possibility that 
other effects may exist whi ch ca n explain the same data. 

Finally, Jaffe et al. claim |T3fij that the notion of the irrelevance of chiral symmetry breaking in the spectrum cannot 
be correct under any circumstances and speculate that hints of parity restorati on in the baryon spectrum are instead 
due to Ua(X) restoration. These observations are at odds with those of Ref. |l28j | and the simple model presented 
here. 

C. D s (2317) and A,(2460) 
1. Experiment 

BaBar's announcement (T2l) of a narrow low mass state seen in the D s n° decay mode has gene rated renewed 
enthusiasm for open charm spectroscopy - in large because of the curious properties of the new state |l3lj . Fig. 
shows fits to the D s ?r mass distributions in the KKtttt and KKitttit modes. The resulting resonance masses are 

M(D S ) = 2316.8 ±0.4 ±3 MeV (155) 



M{D S ) = 2317.6 ± 1.3 MeV 



(156) 
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FIG. 29: Dsn Mass Distribution for (a) D s -> KKty and (b) D s -> KK-kit. 



respectively. The respective widths were determined to be 8.6 ± 0.4 MeV and 8.8 ± 1.1 MeV. Presumably these are 
consistent with detector resolution since BaBar merely conclude that the width of the D s (2317) is less than 10 MeV 
(no confidence interval is reported). 

Evidently the new state is a D s meson which is decaying via an isospin violating mode. Since its mass is below the 
well-known £> s i(2536) and D ,,2 (2573), it suggests that the Z? s (2317) is a J p — + state. This assignment dis allow s 
the decay mode D s (2317) — > D s l, which is consistent with the lack of signal reported in this channel in Ref. jl2l| . 
Similarly, the Z? s (2317) was not seen in the -D s 77 decay mode. 




FIG. 30: The D s (2460) — > D s "/ Helicity Distribution. Solid and dashed lines are Monte Carlo predictions for the J = 1 and 
J — 2 hypotheses respectively^^]. 

A search for the £> s (2317) was also carried out in the D s ir°j mode. No signal at 2317 MeV was observed; although 
a peak at 2460 MeV was com mente d on. The new state was subsequently confirmed by CLEO |l3"2^ in the D*n° mode. 
Further data taking by BaBar [l33j yielded the following improved mass measurements: 

M(L>s(2317)) = 2318.9 ± 0.3 ± 0.9 MeV (157) 

and for the D J2460)[l34j: 



M(D S (2460) D sl ) = 2457.2 ± 1.6 ± 1.3 MeV, 
Mp s (2460) £> s 7r 7) = 2459.1 ± 1.3 ± 1.2 MeV, 



(158) 
(159) 



M(D S (2460) -> D s tt+tt-) = 2460.1 ± 0.3 ± 1.2 MeV. 



(160) 
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The combined result of the latter three measurements is 



M (D s (2460)) = 2459.4 ± 0.3 ± 1.0 MeV (161) 



As seen in Fig. 1301 the D s j helicity distribution strongly prefers a J = 1 assignment for the £^(2460) 139J. Thus 
the low-lying D S (P) quadruplet is complete. 

The following ratios of branching fractions were also obtained: 

1^(2460) ^ s7 ) - 0.375 ±0.054 ±0.057, (162) 



£r(Z) s (2460) -> L> s 7r° 7 ) 

Br(L> s (2460) -> D s n+n-) 
J3r(L> s (2460) L> s tt 7 ) 



0.082 ±0.018 ±0.011, (163) 



Br(D s (2A6Q) ^ D s j) „ „ 

,; < 0.17, (164) 



£r(£> s (2460) -> D s ir° 
Br(D s (2460) -> L> s 7r° 



Br(£> s (2460) -> £> S 7T° 7 ) 

and 



<0.11, (165) 



£?•(£>, (2460) -> L> s (2317) 7 ) 
^(15,(2460) -> D s 7r° 7 ) 



< 0.23. (166) 



The latter three are 95% C.L 

3p 



Recently, Belle |l35l Il36j has reported the observations of the new £> s 's in B decays. The D* (2317) is seen in the 
process B° -> D*^(2317)K _ , D*^j~(2317) -> D+tt°, and the D+^2460) may have been seen (at a much lower rate) in 
a similar chain to D+ 7 . The reported branching fractions are 



Br(B° -> D*+(2317)K") • £r(D*+(2317) -> D+tt°) = (4.4 ± 0.8 ± 0.6 ± 1.1) • 10~ 5 , (167) 

Br(B° -> D*+(2460)K-) • £r(D*+(2460) -> D+ 7 ) = f 0.53 ± 0.20 ± 0.6 +0J6n ) . i(r 5 . (168) 

V —0.15/ 

One may compare these rates to those for B — > D s if with a branching fraction of roughly 4 • 10~ 5 . These processes 
can occur via W exchange or final state interactions. The tree diagram can also contribute to a four-quark component 
of the D s . Assuming that the D s ir° mode dominates the width of the _D S (2317) then implies that the ratio 



Br(B^D s ^K,D)-) ^ 

Br(B a — > D s o(tt, K, D)~) ' V ' 



This is in agreement with expectations for cs states |137j . Alternatively this ratio is expected to be 10 for multiquark 
states. 

Finally, the reported yields for possible charge partners of the D s (2317) are —28 ± 25 for Z?°(2317) and —39 ± 16 
for D+ + (2317). Thus tetr aqua rk models (to be discussed below) of the D s (2317) appear to be ruled out. Also we 
note that angular analysis 138J of Z? s i(2536) — > D* + Ks indicates the presence of a large S-wave amplitude. This 
observation counters the notion that the narrow j q = 3/2 states decay in D-wave. 
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2. Models 

The D s (2317) lies approximately 160 MeV below the majority of model computations. A selection of such model 
predictions and a summary of the experimental masses are given in Table 1X1 VI One sees good agreement between 
data and models for the j q = 3/2 multiplet and poor agreement for the j q = 1/2 multiplet. This is also illustrated in 
Fig. 1311 The poor agreement calls into question the cs interpretatio n of the new D s states; however this interpretation 
has received support from light cone sum rule deca y computations 140] . A comprehensive analysis of the new states 
in the cs model is given by Cahn and Jackson [l4l| . 

The majority of models listed in Table IXIVI predict a D s0 between 2400 and 2500 MeV. The model of Matsuki 
and Morii 148] is unique in predicting the D s q and Dsl masses approximately correctly. The model employs a Foldy- 
Wouthuysen nonrelativistic reduction of an assumed scalar confining interaction and vector Coulomb interaction 13 . 
It is difficult to isolate the cause of the good agreement - it may due to relativistic light quark kinematics or the 
novel assumed spin structure of the static interaction. Further pred ictio ns of this model are discussed in Sect. IIV El 
Finally, Godfrey discusses using radiative transitions to test models [l50| . 



TABLE XIV: Theory and Experimental D s Masses 



D s0 




D' sl 


D s2 


ref 


2480 


2550 


2550 


2590 


GKfl22] 


2380 


2510 


2520 


2520 


ZOPJ142] 


2388 


2521 


2536 


2573 


GJfl43l 


2508 


2569 


2515 


2560 


EOF f 1441 


2455 


2502 


2522 


2586 


LNPJ1451 


2487 


2605 


2535 


2581 


DEfl46l 


2339 


2496 


2487 


2540 


MMfl47] 


2317.3 ±0.4 ±0.8 
2317.2 ±0.5 ±0.9 
2318.5 ± 1.2 ± 1.1 


2458.0 ±1.0 ±1.0 

2456.5 ±1.3 ±1.3 

2463.6 ± 1.7 ± 1.2 






BaBarfl21] 
Belle [139] 
CLEO [132] 



2317.4 ± 0.9 2459.3 ± 1.3 2535.35 ± 0.34 2572.4 ± 1.5 RPP[17] 



Molecules 

Such a strong deviation from expectations point to exotic structure, such as D K molecules, or to a deep misunder- 
standing of heavy-light hadrons. 

Mass (MeV) 

2800 



2600 



2400 



2200 



2000 



1800 

D D* D „ D D' D „ 
s s "s0 S 1 s1 s2 

FIG. 31: The Low Lying D a Spectrum. Lines are theory from Ref.|l8j|. boxes are experiment. The dashed lines show the DK 
and D*K thresholds [2^]. 




The full vector current i s use d. Also, the authors found it was not possible to extract 3 _Di and 3 D2 masses, although recent computations 
have resolved this issue l!49l . 
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As can be seen in Fig. |^ the D s o lies just below DK threshold, while the D s i lies just below DK*. This suggests 
that mixing with thes e con tinua may be important components of the new D s states. Indeed, Barnes, Close, and 
Lipkin have suggested |l5l) that the D s q and D s i may be DK^ molecules. They note that such an assignment 
naturally explains the anomalous masses of the D s (2317) and Z? s (2460). Furthermore, since the D°K + channel is 
roughly 9 MeV higher than the D + K° channel, isospin breaking is induced, thereby explaining the isospin breaking 
discovery mode of the Z? s (2317). The same observation holds for the DK* channel, suggesting that isospin breaking 
may be important in the Z? s (2460) as well. 

The molecular assignment for the D s q can be tested by measuring the radiative decay D s q — > D*j. A simple 
quark model computation gives 1 keV for this El transition 26] . Similarly the D s (2460) transitions are D s \ — > D s ^ = 
7.3 sin 2 (f> keV and D s \ — ► Z?*7 = 4.2 sin 2 (f> keV. Using the heavy quark mixing angle sin 2 </> = 2/3 yields predicted 
rates of 4.5 keV and 2.8 keV respectively. 

Alternatively, the radiative decay of molecular D s j proceeds via the annihilation diagram shown in Fig. 1171 (right 
panel). A simple estimate of these rates gives 26] 



T{D s0 {DK) -> D* sl ) m 35 keV (170) 
T{D sX (DK*) -» D*j) » 45 keV (171) 
T{D sl {DK*) -> D sl ) » 50 keV (172) 

Thus it is possible that such a measurement will serve to distinguish the structure of the anomalous D s states. There 
is also the intriguing possibility of a radiative transition between molecular D s i(DK*) and D s q(DK) states which 
would be driven by the virtual transition K* — > K-y. This rate is estimated to be 17 keV in Ref. [26|, which may be 
compared to the analo gous cs rate of (0.3cos</> + 1.4 sin (f>) 2 1 keV. Recent computations of strong decay rates are 
presented in Ref. |l52fl . 

Finally, if the D s q and D s \ really are molecular states, extra D s j states should appear in the spectrum around 2500 
MeV. Searching for these is of obvious importance. 

Szczepaniak also co nside rs a molecular assignment for the D s (2317) as D s tt bound state with a binding energy of 
approximately 40 MeV |l53l ] . He examines the properties of such a state with an effective lagrangian which incorporates 
the nonlinear realisation of chiral symmetry and Ua(1) symmetry breaking via coupling of pions to the field Q = FF. 
The coupling of D s mesons to pions is affected via a Q 2 D\ term in the lagrangian. Parameters and cutoffs are fit to 
the mass of the D s (2317), and the formalism is used to predict other charged Dgir^ modes, a scalar Dn at 2150-2300 
MeV, and a scalar DK molecule at 2440-2550 MeV. All of these molecules are predicted to have widths of order 
10 MeV. Finally, Szczepaniak also suggests that D*ir and D*K molecules may exist as well. Unfortunately, little 
evidence exists for this collection of predicted molecules. Generally, molecules containing pions are not considered 
viable because the small mass of the pion leads to very small molecular reduced masses, which tends to destabilise 
the system. 

Coupled Channels 

A less extreme version of the molecular picture postulates that the D s (2317) and Z? s (2460) are dominantly cs 
states which are he avily renormalised by mixing with the DK and DK* continua. This approach has been followed 
by Hwang and Kim |l54j who use the Cornell coupled channel model discussed in Sect. II B 31 to compute the mass 
shifts experienced by bare cs states upon coupling to the DK and DK* continua. The authors find that the bare 
scalar D s at 2480 MeV has its mass lowe red t o 2320 MeV, in agreement with the D s (2317). The same philosophy has 
been followed by van Beveren and Rupp p"55l | except that the 3 Pq model is used to drive the transitions between cs 
and the DK channel. Note, however, that both papers do not consider coupling to other channels or renormalisation 
(see Sect. II B 41 for further discussion on this issue). 

Tetraquarks 

In general the D s has a Fock space expansion: 

D s (2317) = z\ ,2 \cs) + Zl /2 \cs{uu + dd)) + Z 1 3 /2 \D°K+) + z\ ,2 \D + K a ) + .... (173) 

The models presented here amount to different estimations of the relative strengths of the different channels. 
Tetraquark models assume Z2 3> Z-\ , Z%. Z4 and impose a, specific structure on the int ernal con figuration. Sev- 
eral groups have pursued this scenario^!! ^i6l |l57l Il58t Il59j ] . Note that this is an old idea |l60lll6l| which builds on 
the proposed light tetraquarks of Jaffe|5ll Il62j . 

A large spectrum of [cq] [sq'] tetraquark states exists and it is natural to attempt identifying some of these with the 
anomalous D s states. The [egjstf] states may be obtained from Eq. 1441 upon making the replacement c — > s. The 
resulting spectrum of Maiani et aL[4!j is shown in Fig. [221 Possible assignments of the Z? s (2317), Z? s (2460), and 
D s (2632) (to be discussed in the following section) are also indicated. 
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Along these lines, Chen and Li|16jj suggest modified four quark configurations of cs(uu + dd) for the Z? s (2317) and 
.D s (2460). Nussinov |l64| suggests that if the D s (2317) is a four quark state then a new isoscalar DD bound state 
with mass less t han 3 660 MeV must exist. 

Browder et aL |l65l | present a modified version of the tetraquark scenario, by assuming that tetraquark states lie 
above 2400 MeV and mix with bare cs states at a similar mass. In this picture, the D s (2317) and D s (2460) are 
predominantly cs states and narrow, while the four-quark states lie above DK threshold and are broad. 

Lattice 

UKQCD 166J have computed the D s spectrum in quenched and Nf = 2 unquenched configurations. Their results 
are shown in Fig. 123 The figure illustrates that it is difficult to come to definitive conclusions on masses in this sector 
due to statistical errors of order 100 MeV and systematic errors due to unquenching of similar size. Notice that the 
D s q mass is approximately 100 MeV higher than experiment - in keeping with naive quark model estimates (see the 
discussion in Sect. II B 4|) . We remark that it is extre mely unlikely that the lattice can differentiate the D s i and D' sl 
states. Similar results are also reported by Bali |l67j . who notes that significant finite heavy quark mass corrections 
to the heavy quark mass limit are present at the charm mass scale, non-degenerate sea quarks may be required, and 
quenched simulations of four-quark states are desirable. 
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3. Summary 

It appears likely that the j q = 1/2 (D s q, D s \) multiplet is difficult to interpret as simple cs states and novel dynamics 
are required. Of the continuum of possibilities, tetraquarks are the most speculative and least constrained by theory 
and experiment. Because of the unknown dynamics, tetraquark models tend to predict extensive spectra. So far there 
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is little evidence for tetraquark supernumerary states and one must regard these models as unestablished. 

Molecular interpretations are similar to tetraquark in the assumption that four quark configurations dominate the 
-D s o and D s \. However the presumed diffuse spatial configurations are very different from those of tetraquark clusters 
(of course, the colour and spin configurations can vary as well) . Sufficiently assiduous experimental study will be able 
to distinguish these scenarios. Typical molecular models also suffer from poor knowledge of the relevant dynamics - 
in this case the nonperturbative gluodynamics that drives meson coupling to the continuum. Thus, for example, if 
S-wave channels dominate the putative DK molecule one expects no charged partners. However, these can exist if 
the DK binding is driven by t-channel exchange. 

The scenario which requires the weakest extension of the naive quark model assumes that mixing of scalar (and axial) 
bare cs states to the DK (and DK*) continua substantially modify the structure of these states. This mechanism is 
thought to be responsible for the anomalous properties of the /o(975) and cio(980) and in general, should be important 
whenever bare quark states couple in S-wave to nearby continua. 1 regard this scenario as the most likely explanation 
of the new D s states. 

Unfortunately, theoretical progress will be difficult due to the previously mentioned unknown gluodynamics which 
drives Fock sector mixing, and the difficulty in consistently renormalising the quark model. But it is becoming clear 
that this effort will be required to understand a slowly expanding subset of the hadronic spectrum. 



D. D s (2630) 

The D s (2630) was discovered by the SELEX collaboration at FNAL(16S| in the final states D°K+ and D s r). The 
peak in the DK mode can be seen in Fig. [21 The combined measured mass is 2632.6 ± f .6 MeV and the state is 
surprisingly narrow with a width of less than 17 MeV at the 90% confidence level. 
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FIG. 34: Observation of the D s (2630) in the DK Mode|l 



The ratio of the partial widths is measured to be 



T{D S -> D a K+) 
T(D S - D a rj) 



0.16 ±0.06. 



(174) 



As pointed out by the SELEX collaboration, this is an unusual result since the DK mode is favoured by phase space. 

This observation has engendered speculation that the D s (2630) is an exotic state. It is unlikely that this state is 
a cs hybrid since the mass of such a hybrid is expected to be roughly 3170 MeV. Possible molecular states include a 
D* s rj system at 2660 MeV or D*oj or D*K* states at 2900 MeV. However the former is a P-wave which is not favoured 
for binding, while the latter are too high in mass to be plausible. 

The large numb er of possible tetraquarks makes it simple to fit the D s (2630) into tetraquark spectra. For example, 
Chen and Li |l37| regard the new D s as a [cs][ss] state - an SU(3) partner to the D s (2317). Claimed destructive 
interference in the decay amplitude leads to a narrow width for the state. The authors predict an axial-vector 
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tetraquark partner of the D s (2630) with a mass M(D S \) w 2770 MeV. Note that the maximally attractive tetra quar k 
channel, which is antisymmetric in col our a nd in spin, is forbidden because of the identical anti-strange quarks 169J. 

Alternatively, Maiani et al. suggest [170] that the Z? s (2630) is a [cd][ds] tetraquark (see Fig. I32|) . They speculate 
that this state does not mix with [cu][us] and hence breaks isospin symmetry. The decay to D°K + would then be 
OZI suppressed. However the state should decay to D s ir° and D + K° with sizable branching fraction. Discovering the 
D s tt mode (along with the observed D s r\ mode) would indicate isospin breaking, as required in this scenario. Maiani 
et al. also predict the existence o f a ch arge +2 state with mass near 2630 MeV, which decays to D s tt + and D + K + . 

Finally, van Beveren and Rupp [l7lT | suggest that the D s (2630) is a heavily renormalised cs state with the aid of a 
K matrix coupled channel model. The model also predicts D* states at 2720, 3030, and 3080 MeV with widths of 2, 
50-60, and 8 MeV respectively. 

If the predicted decay modes or part ner s tates of the D s are not found, the only remaining possibility is that the 
Z? s (2630) is a radially excited cs vector [l72^ (although it is some 100 MeV lighter than the quark model predictions 
of 2730 MeV). The peculiar decay ratio remains to be explained. Experience with the decay modes of the to 
DD, DD* and D*D* points to a possible resolution: transition matrix elements for excited hadrons may have z eroe s 
due to wavefunction nodes. It is possible that such a node is suppressing the DK decay mode. Computation^^], 
shown in Fig. [213 reveals that there is indeed a node but that it occurs at a wavefunction scale which is 20% lower 
than the scale expected in quark models. Furthermore, the DK mode is always larger than the D s r\ mode so that 
explaining the data is difficult in this scenario. 
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FIG. 35: Theoretical Partial Widths for the D s (2630) as a Function of the SHO D s Scale f3 [l 



We have run out of options and must con clud e that th e P .,,. (2630) is an experimental artef act. This conclusion now 
appears likely because searches by FOCUS |l73j |. BaBar p"74[ (see Fig. I36|l . and CLEO p"75| have found no evidence 
for the D s (2630). 



E. P-wave B Mesons 

Preliminary results from CDF and D0 on orbitally excited B mesons were presented at Hadron 05 |l77| . The B 
mesons were reconstructed in J/ipKir and sufficient events were collected to discern the B2 and B\ contributions to 
the rate (see Fig. I57J) . 

The D0 mass results are 

M(Bi) — 5724 ± 4 ± 7 MeV, (175) 



M(B 2 ) - M{Bi) = 23.6 ± 7.7 ± 3.9 MeV. (176) 

Godfrey and Isgurjl^ predict M(B 2 ) = 5800 MeV, 52 ± 12 MeV too high. The Bib are predicted to be "within 40 
MeV of the B 2 ■ Alternatively, the model of Matsuki and Morii [T47]| . which was so noticeably successful for the D s 
system, predicts M(B 2 ) = 5692 MeV and M(B X ) = 5740 or 5679 MeV. 

The total widths, constrained in the fit to be equal (and predicted to be comparable theoretically) are 



r = 23 ± 12 MeV, 



(177) 
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which is consistent with theoretical expectations for the two narrow P-wave B mesons. The experimental resolution 
is estimated to be about 10 MeV. 
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FIG. 37: Qrbit ally-excited P-wave B mesons, as reported by D0. Fitted contributions from Bi — >B-7r, B2 — >B-7r and B| — >B*7r 
are shown 178]. 



F. B c 

The B c pseudoscalar meson was originally seen in the semileptonic decay channels B c — > J/ipivgX at FNAL [l79j . 
Unfortunately the low statistics and unreconstructed final state did not permit an accurate determination of the mass: 



M(B C ) = 6400 ± 390 ± 130 MeV. (178) 



In May 2005 the CDF collaboration 180J announced an accurate determination of the B c mass. The state was seen 
in 14.6 ± 4.6 events in the J/tpn^ mode (see Fig. 1380 with a reported mass of 



M{B C ) = 6285.7 ± 5.3 ± 1.2 MeV. 



(179) 
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Subsequently, the CDF collaboration has also reconstructed 1 1 8 1| the B c in the decay mode B c — ► J/rpl with a 
significance of 5.2 a for £ = fi (or 5.9 a for £ — e). The following ratios of branching fractions have been measured: 



a(B c )Br(B c -> J/^) _ 249 ± 45+ - 107 
a(B)Br(B -> J/^if) ~ u - U45 -o.076 



and 



a{B c )Br(B c -> J/jfe) 
<r(B)Br(B -> J/^iT) 

An updated fit to the £? c — > J/?/;7r decay mode yields [l8ll 



= 0.282 ±0.038 ±0.074. 



(180) 



(181) 



M(B C ) = 6287.0 ± 4.8 ± 1.1 MeV. 
Lastly, the lifetime of the B c has been measured to be 



(182) 



t(B c ) = 0.474. 



I™™ ± 0.033 ps. 



(183) 



This may be compared to theoretical expectations |182j of 0.55 ± 0.15 ps. Remarkably, the inclusive computation of 
this lifetime considers 65 different decay modes. 

The timing of the CDF announcement was arranged with the HPQCD lattice group fl83j to permit them a prediction 
of its mass. The lattice prediction was 6304 ± 12^q 8 MeV. Thus CDF may claim a state, the lattice may claim a 
victory, and the B c appears to carry no surprises. While there is justifiable pride in the present predictions and 
experiment, some pers pec tive may be gained when it is recalled that Godfrey and Isgur predicted a B c of 6270 MeV 
some twenty years ago|l£ 



As stated in Sect. II B 41 this should not astonish, the quark model and lattice gauge theory 
tend to agree on states with a simple Fock space expansion. 



V. CONCLUSIONS 



The moderately heavy charm quark makes charming spectroscopy especially interesting since it lends confidence 
that the behaviour of soft QCD can be extracted from the data, yet it is not so heavy that one is simply exploring 
an analogue of the positronium spectrum. Thus the new charming mesons have generated much enthusiasm. The 
A(3872), D,,(2317), and D s (2460) remain enigmatic and highlight a less than perfect understanding of QCD. Can we 
be seeing the expected breakdown of the naive quark model? Where and how do gluonic degrees of freedom manifest 
themselves? Have we entered a new era of 'mesonic nuclear physics' ? 

If the V(3872) is not canonical charmonium then the simplest explanation lies in expanding the Fock space de- 
scription of the state. With 1 ++ quantum numbers this is likely to mean ccqq. Whether a diffuse four-quark state 
dominates the X, as in the molecular interpretation, or a compact state dominates as with tetraquarks, or the four- 
quark component heavily renormalises an erstwhile x'ci ^ s a matter of dynamics. This dynamics can be subtle, but 
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often leads to different and testable predictions. Thus, for example, an entire spectrum of tetraquarks is predicted 
while a modest number of additional molecular states are expected. Alternatively, the idea that the X is a purely 
kinematical cusp effect can be checked by measuring the rate for B + — > K°D + D*°. If no enhancement is seen in 
this channel it implies that significant final state interactions are creating the X signal. Finally pinning down the 
DDn mode will be vital in establishing the magnitude of cc-ccqq mixing. And this will tell us something about the 
gluodynamics which generates mixing. 

The X, Y, Z charmonium complex near 3940 MeV poses interesting problems. Prima facie, it has some of the 
expected features of a 2P multiplet. However, if the Z is interpreted as the x' C 2 an d the A(3940) is the x'ci then the 
multiplet displays an inverted spin orbit splitting. This unusual situation would then highlight significant confusion in 
our understanding of the Dirac structure of confinement. Alternatively, the unusual u>J/tp decay mode of the F(3940) 
has led to suggestions that it is a charmonium hybrid. This suggestion does fit with (rather shaky) expectations for 
hybrid decays; however, the mass is much too low compared to (somewhat less shaky) quenched lattice computations 
and prudence - and more data - are called for. The evidence for the Y(3940) is not especially robust and one must 
hope for additional experimental effort in determining the properties of this state. Lastly, if the X(3940) is the x'ci 
then why is the Xci not apparent in the same data set? 

The K(4260) is an intriguing candidate for a vector charmonium hybrid. This possibility may be pursued by 
searching for other S + P decay modes such as XcO^ and for nearby pseudoscalar and J PC = 1 h hybrids. The 
resonance interpretation of the Y may be confounded by the proximity of the D\D channel which will cause an 
enhancement due to the threshold cusp effect and possible strong D\D rescattering. This possibility can be tested by 
measuring analogous strange systems such as D 8 \D and D S \D S . Finding manifestations of glue in the spectrum is a 
primary goal of nonperturbative QCD and one must hope that the Y will be studied in greater detail. CLEO-c may 
be able to assist substantially in this effort. 

The masses of the D s o(23I7) and D s i(2460) point to a significant lack of understanding of QCD dynamics in the 
cs system. Possible interpretations include tetraquark states or molecular DK or DK* states. The scenario which 
requires the weakest extension of the naive quark model assumes that mixing of scalar (and axial) bare cs states to the 
DK (and DK*) continua substantially modify the structure of these states. An alternative explanation is provided 
by the chiral doublet model, but it has been argued that this model is likely not relevant to low-lying heavy-light 
states in Section IIV B 21 Firmly establishing the analogous D st ates will assist in understanding the novel dynamics 
which appears to be present in P-wave heavy- light hadrons |l84) . 

The advent of the pentaquark states and the new mesons has led to increased interest in diquark correlations 
as effective degrees of freedom. However, the subsequent collapse of the pentaquark states raises serious concerns 
about the reliability of chiral soliton and other pentaquark models. Similarly, an infinite collection of tetraquark 
states has been predicted. However, there is no evidence for any of the predicted low lying states. Clearly a degree of 
circumspection is called for. 

Alternatively, many of the proffered explanations of the unusual states involve mixing with higher Fock components. 
These explanations are economical in that new dynamics or degrees of freedom need not be postulated. However, they 
must be carefully constructed because it is relatively easy to destroy previous agreement with experiment. A simple 
organising principle seems to be emerging: the presence of nearby S-wave meson-meson thresholds can strongly affect 
meson structure. While this observation can serve as a useful rule of thumb, turning it into a quantitative theory is a 
daunting task because it is driven by unknown nonperturbative gluodynamics. Of course, understanding gluodynamics 
is one of the great remaining challenges in th e development of the Standard Model, and the effort is worthwhile. 

More than a decade ago Isgur stated |l85j . "Twenty five years after the birth of the quark model, light quark ... 
spectroscopy is in a deplorable state" . While the situation is slowly improving, it is the formerly comfortable heavy 
quark sector which is now under siege. And new data seem to reveal new ignorance with each passing month. It 
is clear that the simple constituent quark model must fail somewhere - are we seeing this? The new heavy meson 
spectroscopy challenges our understanding of QCD. Can we meet the challenge? 
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APPENDIX A: BOUND STATE DECAY FORMALISM 



While the use of the weak binding decay relations of Sect. Ill D II are undoubtedly accurate for the X(3872), it 
is useful to examine molecular decays in more depth to assess the accuracy of the weak binding approximation. 
Furthermore, it was assumed that the p J/ ip decay to irnJ/ ip proceeded as if the p J/ ?/> channel belonged in the bound 
state instead of in the continuum. This simplifying assumption will also be justified in the following. 

We first consider the generic strong decay of a molecule, X to a final state a(3B via a bound state component, AB 
illustrated in Fig. [23 In the 3 P model this amplitude is 
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FIG. 39: Generic X -> (AB) -> af3B Decay. 



A(X^AB-A^al3) = Z^ x(AB) ( qa + q0 ) J A (fc)^(fc + | - ^ • (k- ^ + | ) (Al) 

where (tr) is the matrix element of the Pauli matrices in the Aa(3 spin states and q a and qp are the momenta of the 
mesons a and /3 respectively. 

The differential decay rate is given by 

dT = \A\ 2 dfo = Z AB \cj> x{AB) {q a + qp)\ 2 A{A a/3)| 2 (A2) 

where d03 is the three particle phase space. Notice that the use of the 3 Po model is no longer relevant. Integrating 
and using the recursion property of phase space (Eq. 37.12 of Ref. 01) gi yes 

T(X -> a 0B) = Z AB J jty-\0 x{AB) ( q )fr(A -> a/?) (A3) 

where the A meson is taken to have energy E A — M(X) — E B (q) and momentum q. In the weak binding limit the X 
wavefunction is strongly peaked at the origin and one obtains 



T(X -> a(3B) w Z j4B r( J 4 a/3) 



(A4) 



where now T(A — > a/3) is the experimentally determined free space partial width of the A. This identification requires 
M(X) — M{B) w M(A) which is of course also true in the weak binding limit. Thus the weak binding formulas are 
justified for the case where the channel AB is a bound state component of the molecule. 

The use of the weak binding decay relation for X — > ittt J/tp is more problematic because the pJ/ip channel should 
be considered part of the continuum rather than part of the X wavefunction as assumed in Sect. Ill D II Indeed, this 
process is second order in perturbation theory, as illustrated in Fig. 1401 rather than first order as above. In this figure, 
the boxes represent intermediate states which are summed. 




FIG. 40: Second Order Transition X(AB) -> CD -> a/3D. 
The perturbative expression for the decay amplitude is 



A(X(AB) -^CD-> apD) = T(q a + q p ) 



A{B -> a(3)) 
M{X) - E CD + 



(A5) 



where 
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T(q) = J ^^(I)x(ab){p)Tab^cd{p -> ?)• (A6) 

and T is the perturbative scattering matrix representing the subprocess AB — > CD in Fig. 1401 
Integrating over three body phase space yields the following expression for the width 

r <* - °*°> = / H r ' c - ■*> mm - ecmw + r^4 < A7 > 

where the total width of the C meson has been restored to the propagator and it assumed that the D meson is stable. 
The width for the C is to be evaluated at the kinematical point Ec = M(X) — Ed(Q). 
In the limit where the C particle is narrow this approaches 



) E £ E_ 

' M(X 

as expected. Furthermore, if the secondary amplitude C — > a/3 is independent of momentum Eq. IA7I simplifies to 



T(X - CD) = 2nQ^^Br / \f \ 2 dtl (A8) 



T(X — » CI? — > a^D) = T(C - a/3) | — — - = Z CD r(C -> a/3). (A9) 

Thus the weak binding expression is recovered under the conditions that the production side of Fig. 1401 and the decay 
side factorise. In general one expects no dependence of the C width on the momentum of the C so the constraint is 
certainly plausible. 

APPENDIX B: PION EXCHANGE POTENTIALS 

The standard quark-pion interaction (as used in nuclear physics) is given by 

— ' <2 3 :e^7 m 7 5 tV> • d^n (Bl) 



V2U 

where f„ — 92 MeV and g is the pion-quark coupling constantly. Effective pion exchange potentials are derived 
from the nonrelativistic limit of the appropriate t-channel scattering amplitude where the vertices arc computed from 
the pion-quark interaction given above. 

Extraction of the nHH vertex is simplified by using the Hamiltonian |l86| 

where the sum refers to a sum over quark flavours. The computations may be simplified further by considering forward 
propagating neutral pions: 

For example, the effective one-pion-exchange induced NN interaction is obtained from the ttNN vertex which is 
given by 

(p^H\ P] -) = ^-l-kn . (B4) 
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TABLE XV: One-pion-exchange NN Interactions 



isospin J 



channels 



potential 



1 

1 0" 
1" 

1+ 3 



1+ 3 Si/ 3 D 1 (the deuteron) Eg. IbTTI 



-fV C(r) 
f V {C(r) - 4T(r)) 
25VoC(r) 



This relation has been derived using the proton wavefunction 



\Pt) = 7/g (K w i d t> + NMt) - 2|w T u T d x )) . 



(B5) 



The NN scattering amplitude is then given by 



25 / g \ _ _ a i ■ ka 2 ■ k 



9 V72/ w 

which yields the relationship gNN-n — § 9m ^^ . 

Fourier transforming gives the generic pion-induced interaction 



k 2 + ml 



(B6) 



£,V fi ■ t% [<7i • a 2 C(r) + (3(Ti • r<?2 ■ f — o"i • (T 2 )r(r)] 



(B7) 



where £ is a numerical factor determined by the channel in question and 



V 



247T/2 



1.3 MeV, 



(B8) 



C(r) 



(B9) 



and 



Tfr) 



+ 1 C(r). 



(BIO) 



Nucleon-nucleon interactions may now be evaluated using £ = and (ri • f^) 
TViV interactions are as listed in Table IXVl 



21(1 + 1) - 3. Thus all possible 



V 



25 



K 



C(r) 





n/8 



V8 

-2 



T(r) 



(Bll) 



This corrects errors in Table 1 and Eq. 30 of Ref . |6j| . 

The interaction is often regulated in an attempt to mimic unknown short range dynamics. It is traditional to do so 
with a dipole form factor such as (A 2 — m\ )/ (A 2 + 1) . The resulting softening of the potential is crucial to reproducing 
deuteron properties. In principle, the form factor can be computed by convoluting the spatial pion and nucleon 
wavefunctions. However, attempts to do so yield unsatisfactory results and the dipole form factor has been retained 
in X the molecular computations reported here. Finally, we note that the form factor induces a very short range 
repulsive core (figures 1 and 2 of Ref. [6l| are misleading) that, in part, explains the phenomenological usefulness of 
this method. The cutoff has the further benefit of providing a factorisation scale between long range pion exchange 
and short range quark exchange dynamics. 

Computations for meson-meson interactions proceed in a similar way. For example vertices for K or K* scattering 
are given by 



G2 



and 



TABLE XVI: One-pion-exchange MM Interactions 



state 


isospin 


J 


channels 


potential 


DD* 





n — h 


3 n 

Po 


-3Vb(C(r) + 2T(r)) 


DD* 





0~~ 


3 7"> 

Po 


+3Vo(C(r) + 2T(r)) 


DD* 





1 + 


-Pi 


- Wo (C(r) -T(r)) 


DD* 





1 


3 7"> 

Pi 


+3V (C[r) -T(r)) 


UU 


u 


1 






7—1 rS* 
DD 





i H — 
1 


3 c* 3 n 


—o- (Eq|B14J 


DD* 


1 


n — h 


3 D 


+Vo(C(r) + 2T(r)) 


DD* 


1 


0~~ 


3 D 

-Po 


-Vo(C(r) + 2T(r)) 


DD* 


1 


1 h 

1 + 


3 7"> 

Pi 


+Vo(C(r)-T(r)) 


DD* 


1 


1~~ 


3 n 
Pi 


-V {C(r) - T(r)) 


DD* 





1 ++ 


3 ft, 3 Di 


— fEalBTH 


DD* 





1 A — 


3 ft, 3 -Di 


+ (Eq|B14J 


pi* i—,* 

L> D 





u 


OO, -L<0 


■ (EalB15[l 


D*D* 







3 n 

Po 


-3Vb(C(r) + 2T(r)) 


D*D* 





1 ++ 


Di 


+3l / o(C(r) - T(r)) 


D*D* 





1 — 4- 

1 ^ 


3 7-i 

Pi 


-3Vo(C(r) - T(r)) 


D*D* 







3 ft, 3 -Di 


3- CEalB14^ 


D*D* 





1" 


1 Pi, 6 Pi, 5 Fi 


6- (EqlElll 


D*D* 







1 D2, 5 ft, 5 -D2, 5 G*2 


6- fEalB17^ 


D*D* 





nH 


3 7-» 

D 2 


-3Vo(C(r) - T(r)) 


D*D* 


1 


n++ 


L So, j Dq 


-2 ■ (ha IB15H 


D*D* 


1 


o-+ 


3 Po 


+Vo(C*(r) + 2T(r)) 


D*D* 


1 


1++ 


5 d 


-V (C(r)-T(r)) 


D*D* 


1 


1-+ 


3 Pi 


+V (C(r)-T(r)) 


D*D* 


1 


1+- 


3 ft, 3 P»i 


-(Eqinnj 


D*D* 


1 


1~~ 


1 Pi, 5 Pi, 5 f\ 


-2- (EalBT6l 


D*D* 


1 


2++ 


1 £ ) 2, ft, P>2, G2 


-2- (EqlEJI 


D*D* 


1 


2+- 


3 D 2 


V (C(r)-T(r)) 



(K\H\K) = -jLj-k^ 



(B12) 



(K\H\K) 



9 k^. 



(B13) 



Forming the t-channel scattering amplitude and Fourier transforming then yields the effective one-pion-exchange 
meson-meson interactions shown in Table IXVll 



V = -V Q 



1 
1 



C(r) 





-V2 



-V2 
1 



T(r) 



(B14) 



r 



-To 




(B15) 



V 



-Vb 



1 
-i 






-I. 



C(r) 



( 



4 

7 
5 

2.-, 



J5_ 

25 

■ S / 
5 / 



T(r) 



(B16) 
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